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ABSTRACT 


In the present work, two-dimensional laminar com- 
pressihle boundary layer problem under the assumption of 
either Prandtl number equal to unity or that of low mach 
number, has been solved numerically using Parametric Diffe- 
rentiation Method. At the end, a general program based on 
Fortran IV has been developed by which through the use of 
the initial value obtained by the Parametric Differentiation 
Method as starting guess, the original non-linear equation 
can be solved to any degree of accuracy. 
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CHAPTER 1 


IHTRODUCTIOH 

The equations governing the two dimensional compres- 
sible laminar boundary layer are coupled in addition to being 
non-linear. Except for the special case of a flow with Prandtl 
number unity over an insulated surface, these equations are 
difficult to solve. 

With no standard method can the solutions of these 
two point boundary value problems be expressed in a closed 
form. The following two methods are generally resorted to 
for solving these equations. 

1. Forward integration 

2, Integration by successive approximation 

Both of these methods are highly iterative and a good 
amount of guess work has to be done for rapid convergence to 
the solution. The aim of the present work is to dispense 
with the guess work and to solve the equations by parametric 
differentiation (Ref. 5). As we shall see, this method reduces 
the equation to solved to a linear system which, though of 
higher order, can be solved by a predetermined number of 
iterations. 

Satisfactory results obtained in the case of the com- 
pressible boundary layer encouraged us to develop^ a general 
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program which can handle all those ordinary differential 
equations which can he solved by parametric differentiation 
method. ' In the last chapter we have presented this general 
program. Due to the flexibilities incorporated in the pro- 
gram, it can handle, in addition to the reduced linear sys- 
tem of equation, the parent non-linear system of equations. 
This flexibility is used when a high' accuracy is needed. Here 
the solution is obtained in two steps in the first step, the 
usual solution is obtained using parametric differentiation 
method. In the second step, the nonlinear equations are 
directly solved. Here the missing initial boundary conditi- 
ons are taken from the solutions obtained in the first step. 
Is initial boundary condition taken are very close to the two 
one. This method converges very rapidly in a few iterations 
to the ti.ue solution . 



CHAPTER 2 




POEMULATIOH OP THE PROBLEM (Ref. 4) 


2.1 Stewart son’ s Equation: 

steady two-dimensional compressible laminar boundaiy 
layer for .perfect g/ ' "• are given by 


"5F ( + T7 




9u 

J 




( f v) = 0 (Cont. Eqn. ) 

9p 


( 1 ) 

( 2 ) 


U ^ + -1. ( A 

^ dx 8y ^ Pr ay 


) - 


3r 

(Energy Eqn. ) (3) 


where , 

p = pressure 
P = density 

u = longitudinal velocity 
V = transverse velocity 
X = longitudinal coordinate 
y = transverse coordinate 
h , = enthalpy 

^ = coefficient of viscosity 

Pr = Prandtl Humber 

The viscosity law to be assumed is 



( 4 ) 
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The constant ^ is used to match the viscosity with the 
Sutherland value (Ref. -.3) at a desired station. If "this station 
is taken to be the surface assumed to be at constant temperature 
the result is 

vrrTFT (^044^) ^ 5 ) 

w su . 

where 

t^ = wall temperature 

t^ = free strain stagnation temperature 
~ Sutherland's constant 

2.2 Stewart son' s Transformation: 

Velocities in the equation of motion (1) to (3) can be 
expressed in terms of derivatives of stream function as 


fy = ( PU/, t P 


Now introducing transformation 

a 


X 


X P^ 

r 

p 

o -^0 


dx 


• 3s 
■nr 


_e 

ro 


dy 



( 7 ) ■ • 


(Subscript e refers to local condition at the outer edge of 
the boundary layer (external). The subscript o refers to the 
free-stream stagnation values.) Where a is sonic velocity- 





Me get the following equations: 


+ Vy = 0 


UU 


X 


+ Si„ a 


e eX 
S, 


o "YY 


USj+ vs^ = (-g 


1-Pr 
Pr 1 


(( Y -1 )/2)l 


,^2 ■. 
e / 


l 1 + ((Y -1)/2)M 


(.JL)' 


YT 


( 8 ) 


(9) 


whe re 

y = ratio of specific heats 

S = h/f * 

S = (10) 

a is local sonic velocity. 

s is enthalpy function 

and h is local stagnation enthalpy, 
s 

Transformed velocities U, V are related to stream func 
tion hy the following equations; 

u = -% , y = 

The transformed longitudinal velocity U is related to the 
longitudinal velocity in physical plane hy 
IJ = (a^/a^) u 

The boundary conditions applicable to eqn. (7) to (9) are 
U(X, 0) = 0 

V(X, 0) = 0 

s(X, 0) = s^ or ^ (X, 0) = (|f)^ ('1'') 

lim s = 0 

Y-oo 
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Now under the transformation 


■f = f('T) ) rT2 4oUgX)/(m+1) 

'i = Y ro:5 (m+irUg7Ti^xT 


( 12 ) 


The above system of partial equations is reduced to following 
ordinary diffrential equation 


f + ff = p(f ^ - S) 


S + Pr f S = (1-Pr) 




I 1 + ( (^-1 )/2)M^ 


, t Iff HO 

(f f + f (1 


The pressure gradient parameter |3 is defined as 


where m is given by 


2m 

m+1 


where C is’ constant. 

The velocity ratio U/U. = u/u = f* 

© © 

In above equation prime denotes differentiation with 
respect to . 

The boundary conditions are 

f(0) = f’(0) = 0 

S(0) = S 

f 

lim f = 1 


( 14 ) 
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Now in tbe energy equation, the right handside is not 
functionally consistent with arbitrary Mach no. Me which is 
a, function of X. 'Hence to be consistent with left handside 
right handside - of energy equation must be either constant 
or function ofY| . Now this is true for following cases 

1 . 'Then external Mach number is constant as is in the 
case of flate plate, 

2. , External Mach number is negligibly small. 

3. Prandtl number is equal to unity 

4. Specific heat ratio is equal to unity. In reality i-.. 
cannot be equal to uhity. 

5. Mach number is very high (Me -» «>) as in this case 
right hand side will be equal to 2. 

Here problem will be treated assuming either Prandtl 
number is equal to 1 or Mach number is very small. 

Hence system of equations to be solved are 

f + ff = p (f ^ - s) (15) 

s + fs Pr ~ 0 ( 16 ) 

with boundary conditions (14) 



CHAPTER 3 


]yiETHOD OE SOLUTIOI 

3. 1 Reduction of Nonlinear Equations to linear Porm: 

Essence of method based on parametric differentiation 
'is to differentiate the original equation ■#itb respect to a 
parameter, and to solve the resulting linear equation. 

Bow to apply this method we have to know the solution 
of the equation for one value of parameter and then we march 
forward for other values of the parameter. 

Now differentiating equations (15) and (16) with res- 
pect to p we get, 

11 f ff 11.19 .11 . . 

G + Gf + fG = (f ^-S) + p(2f G -T) (17) 

t" + Pr ?«.(fT' + GS') = 0 (18) 

with boundary conditions 

g’( 0) = G(0) = 0, T(0) = 0,T(°o) = 0 

where 

df/d6 = G 

(19) 

dS/dp = T 

wHere prime denotes differentiation with respect to 

Equations (17) and (18) are linear in G and T. 

Now after solvinfc for f and S for initial value of 
p, we substitute, for f and S in linear equation (17) and (18) 
and evaluate G and T. Prom these values of G and T, we 
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solve for f and a from equations (19). This cycle is repeated 
till p reaches the desired value. Both the systems of equa- 
tions (17, 18 ) and ( 19 ) are solved usin^r Runge-Kutt a method 
with proper increment in y and p respectively. 

Initially such a value of parameter is chosen that at 
this value either the solution is known or it can be obtained 
very easily. 

Here the initial value of parameter p is chosen as zero. 
For this value of p, the equations (15 , 16) are reduced to 

f + ff = 0 (20) 

s" + Prfs' = 0 (21) 

with boundary conditions (14). 

Equation (20) is Blasius equation whose solution is 
known .or can be obtained very easily. 

Solution for other values of p .from that at the initial 
value of p is obtained following the method outlined above. 

In actual numerical solution, increment ini] has been 
taken equal to .05 and that in p. equal to 0. 1 for favourable 
pressure gradient and -0.05 for adverse pressure gradient. 
Solutions were obtained for various Prandtl numbers and wall tern- 
pera.tures at different values of p. 

We have seep equations, 'to. -be solved for p equal to zero 
are uncoupled and, simple but still it is nonlinear. We can 

“Pil T’+'hoT' .cs TTTi'nl 1 "Flir “hVio om ’Htt noo rt-P *naT»amp"f:T'i O, 
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differentiation method. In Blasins equation there is no 
parameter hut we can artificially introduce a parameter as 
shown helow. 

let us consider the equation, 

f’” + (1-P)f" + P ff" = 0 (22) 

with houndarj'- condition 

f(0) = f'(0) = 0 

f (°°) = 1 

Above equation for P equal to one is reduced to Blasius 
equation, Bor P equal to zero equation is reduced to, 

TH n 

f + f = 0 (23) 

with boundary condition 

■ f(0) = f' (0) = 0 f’ (~) = 1. 

Solution of equation (23) is very simple and is given by 

f ~ ^ 

Now differentiating equation (22) with respect to P, 
we get, 

w"'+ (1-P)w”, + f" + ff" + PWf" + Pfw" = 0 .(24) 

with boundary conditions 

W(0) = w' (0) = 0 

W(°°) = . 1 


( 25 ) 
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Ttoere, 

df/dP = W (26) 

Equations (22, 24 and 26) are solved in the same way as 
equations (17, 18 and 19). 

Equation (23) has been solved with increment in 
equal to 0.05 and that in P equal to 0.05. Besults obtaiiied 
tally fairly with the standard result. Results are shov/n in 
table (1) oa the next' page. 

Similarly initial solution for both the equations 
(15, 16) can be obtained by assurainy suitable equations. 

This has been illustrated in (^Chapter 5. 
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Table 1 ' 

Solution of Equation (22) by Parametric Differentiation 
for Pr = 1r (Numerical solution for Blasius Equation) 


ETAfV]} 

P 

ETA (f j) 


0 

0 

3.8 

2-583 

0.2 

.0094 

4.0 

2.782 

0.4 

.0375 

4.2 

2.982 

0.6 

.0843 

4.4 

3,181 

00 

• 

o 

.1495 

4.6 

3,381 

1.0 

.2328 

4.8 

3.581 

1.2 

.3333 

5.0 

3.781 

1.4 

-4503 

5.2 

3.981 

1.6 

.5824 

5.4 

4.181 

1.8 

. 7282 

5.6 

4.380 

2,0 

.8860 

5.8 

4.580 

2*2 

1.054 

6 . 0 

4 . 780 

2.4 

1.231 

6.2 

4.980 

2.6 

1.414 

6.4 

5.180 

2.8 

1-602 

6 ,6 

5.380 

3.0 

1.794 

6.8 

5.580 

3.2 

1.989 

7.0 

5.780 

3.4 

2. 186 

7.2 

5,980 

3.6 

2.384 

7.4 

6.178 



CHAPTER 4 


RESULTS OBTAINED 

4.1 Velocity and 'Enthalpy Function; 

The velocity and enthalpy functions are presented in 
tabular form. 


The distance y normal to the surface in the physical 

plane is related to the similarity variable ir| through 

equation (6) and (12) as (Ref. 4). 

p ■ a 
12 . _i 
Pc ' S' 


y 


I’O 


e -4 m+ 1 


TT ^ t '■ 

e o 0 


where 


t/t^ = (1 + 


‘tJ-I 


if ) S f'^ 

e ^ 2 e 


(27) 


2 " e ^ ^ - 2 

4. 2lntegral Thicknesses: 

The boundary layer integral thicknesses in the 
transformed plane are defined by the following relations (Ref. 4 

Displacement thickness: 


tr 


X 


9 


tr 

X 


Thermal thickness; 


m+1 

U X 
e 

4 -T- 

4 

O 1 

ckness: 


i m+1 


4 2 

Vq 


/ (S - f ) 
0 


(28) 


f f (1 - f ) d 
o 


(29) 


.i sil 

2 ;4 ^io' 


f S dr| 


(30) 
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Convection thickness: 



m+l 

2 ” 


N O 


/ (s-1) f’a = -s; 

O 


(31) 


Numerical values of these thicknesses s-re evaluated for 

different values of 6. Pr and S . 

w 

4,3 Shear and skin friction: 


The quantity that is of primary interest in boundary 
layer calculation is the shear stress at wall which in 
nondimensional form can be given as (Ref, 4) 



1 

2 


fw "e 



4'. (2 .'(1+sp) 


\ mjj ^ 0 

i "2" V 



which can be rewritten as 


'C „ iT^ 


w _ 


w 


m+1 d I n X 
2 d In ’ X 


(33) 


where is called local skin friction coefficient. 


Re 




w 




w 


w subscript denotes properties at wall. 


4.4 Heat Transfer: 

S* = ds/d'jo at wall corresponds to heat, transfer 
across boundary layer. This is related to stagnation enthalpy 
derivative in the physical plane by 

VlD,+.1 


h. 


( «) 


( 


f ^ 


.l^nrrnm 


) S 


(34) 
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A non-dimensional quantity in connection with heat 
transfer can he introduced (Ref. 4) as 


Tfu = 


x(8t/6y) 


w 


^ ~w ^ iro” 1 Ei+1 d In X 
- V 4 “2" d Tn x 


t - t 
0 w 


(55 


Reynolds analogy; 


A simple modified Reynold analogy can he defined as 




(-S, 


/“s.' 


(3S 


w 


• 1 ) 


This quantity is the reciprocal of usual Reynold analogy 
quantity. This factor is tabulated for different values of 
p at different yalues of Pr and Sw , 

Numerical results are given in Appendix 1. Comnuter 
program for this is given in Appendix 2. 


4.4 Some Comments on the Result: 

In Pig. (1, 2, 3)velocity profiles have heen shown. We 
observe that 

i) for a given wall temperature and Prandtl numher> the 
initial slope of velocity increases with the pressure gra- 
dient parameter (Pig. 7 also). 

ii) for adverse pressure gradient, an inflexion occurs 
within the boundary layer that moves outwards with increase 
in magnitude of negative pressure gradient parameter. 
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iii) for Sw > 1 , (or for heated surface) under favourable 
pressure gradient a velocity overshoot (Fig. 2) occurs within 
the boundary layer. This overshoot increases as the pressure 
gradient parameter hecames more favourable.. The over shoot 
takes place due to the fact that when wall is heated in a 
favourable pressure gradient, the density in certain region 
of boundary layer is lowered so that in spite of viscous drag, 
the flow there is accelerated more than the external flow by 
the pressure force. 

iv) velocity boundary layer increases with increase in 
Prandtl number. 

v) In Pig. (3 , 4 , 5)otagnation temperature profiles have 
been drawn. We see ' for a given Prandtl no, and specified 
wall temperature profile vary very slowly for different pres- 
sure gradient parameter£. This variation becomes prominent 
with increase in the wall temperature. Like initial velocity 
gradient , slope of initial temperature gradient at wall in- 
creases as the pressure gradient becomes more favourable 
(Pig. 8 also). 

vi) Thermal boundary layer thickness decreases with in- 
crease in Prandtl number. 

By inspecting tabulated results, we obterve that 
for low wall temperature in favourable pressure gradient 
displacement thickness is negative. This is due to the fact 
that the fluid in contact with the cold wall has higher den- 
' sity than that in the external flow. Hence more fluid /area. 



CHAPTER 5 


SOME lETAILS OP COIVIPUTER PRO&EAMTgNG 

5-1 Iteration Scheme Used for Satisfying Boundary Conditions 

Before explaining about the iteration scheme, we will 
mention one of the important properties of the ordinary linear 
differential equation. 

Let the system of linear differential equation be given 
by 

L(D) = 0 

How this can be reduced to a set of linear differential equa- 
tion of order one, number of equations of order one being equal 
to the sutri of the ordeisof the original system— equations. 
Any derivative of any dependent variable in the original equa- 
tions will correspond to one dependent variable in this new 
set cf equations of order one. 

Let x.| , Xg, ... be dependent variables in the new 

set of equations. Let x^(0) , x^(0),... ^p(O) missing ini- 

tial conditions whose value we have to know for solving the 
problem as initial value problem. , 

Let the conditions to be satisfied at the end point be 

Xg_(R) = EM(1) 

■x^(R) = EH(2) 
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where M(l), M( 2)5 ... M(n) are values of variables x , x etc, 

3. c 

to be satisfied at the end point R. Let us relabel the missing 
initial conditions as &(1), G-(2), ..,G(n). 

Thus 

G(1) = x^(0) 

G(2) = X;^(0) 


G(n) = Xp(0) 

Similarly, x , x , x are ^felebclled as 

cX> O ,j 

Y(1) = x^(R) - M(1) 
y(2) = x^(R) = RN(2) 

Yln)= x^(R) = M(n) 

let G(i) = C(i), for i = 1, n 

where C(i) is a value assumed for G(i). 

Row for linear equation it is well known that if all the 
initial conditions are kept fixed, and an increment is given to 
the initial value of a variable then 

^Yi(3) = 0.^ had) (57) 

where Ay^(qj) is increment in 1 ( 2 ) due to the increment G(i) 
in G(i), all other initial values, of variables being kept constant 
6. . is a constant. 

Graphically it is shown in the adjacent figure 
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Now if there are increments , x^ in &( 1)5 G-(2) G-Cn) 

then total increment in Y(o) will he given by 

Ay(3) = 2 ay (o) = S C X. (38) 

i=1 ^ i=1 ^ 

Now we can treat the problem inversely in the following manner. 
Giving different values to initial conditions, we can calculate 
coefficients like C.. from (37), 

Now if AY(j) is departure of Y(,'i) from, the prescribed NN(3), 
when initial values are assumed as 

G( i) = C( i) , i = 1 , n 

Now for a given AY(l), we can solve for x^^ from the 
simultaneous algebraic equations (38). After obtaining x^, correct 
initial conditions which will satisfy the given boundary condition 
will be given by . 

G(i) = C(i) +• x(i)-,- i = 1, n (39) 

Now above is strictly true for linear equations but 
approximately true’ for non-linear equations where guessed initial 
values are near the true values. 

In above method of interpolation for obtaining correct 
initial values, effect of modifications of all initial values 
on any particular final value has been taken into consideration. 
Thus in non-linear case it , can be hoped that convergence to the 
correct' initial value will be obtained more rapidly in the scheme 
of , interpolation than that can be obtained in a scheme where into 
polation is done only between one initial value and corresponding 
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5 . 2 Example 

General program can best be imderstood with reference to 
some equations to be solved. Let us take the equations (15) 
and (16) as an example. 


( 40 ) 


Consider the following equation 
f"’ + B (ff' + (s - f'2) p) + f” (1-B) = 0 
s" + B fs' Pr + s' - Br' = 0 
with boundary conditions f ' (O) = f(0) = 0^ f ’ (°°)= 1 

s (0) = Sw, s(°°) = 1 

Above set of equations is reducible to equations (15? 16) 
for B s^equal to unity. Now for B = 0, equation (40) is reduced 

to a very simple form (41) 

fttf + f»'(i_B) = 0 

s " + s ' = 0 

whose solution is given by 

f = 1"! + e~^^ -1 ■ 

s = 1 + (Sw - 1 ) e" 

Erom the solution of equations (40) for B equal to 9^zero? we 


(41) 


(41a) 


can ma.rch forward by parametric, differentiation method to get 
the solution at B. equal to unity. 

Though equation (40) looks more comnlicated than the 
parent equation nonetheless it is chosen as it gives the startin 
solution readily. 

: low differentiating (40) with respect to B we get 



21 


G"'+(ff'+(s-f'^) +B(Gf»+fG"+(T-2f' )p)-f”-G"( 1-B) =0 


BT’-s’ = 0 


(42) 


(43) 


T"+ fs' Pr + B Gs' Pr + BfT'Pr + T' - 
with boundary conditions 

G(0) = G'(0) = G'(co) = 0 
and T(0) = T(°°) = 0 

where 

r = M m - bs 

^ dB ’ ~ dB 

Parent equ8.tions are 

f'"+ ff" + (s-f’^)p = 0 

s '' + fs ' Pr = 0 

with boundary conditions 

f(0) = f'(0) = 0 
f '(«=)= 1 

and s(0) = Sw, s(°o) = 1 

Now we have 3 sets of equations. 1st set is equation (41), 

second set is (42) and third one is eqn. (43). In PDM (Parametric 

Differentiation Method), we first solve 1st set, then second set. 
If higher accuracy is required we solve the third set of equations 
directly. In solving the third set we obtained , starting values 
by solving the second set of equations by PDM. 

In this program any dependent variable is represented as 

■ f 

double subscripted variable ,Q(n, m) where n denotes value of 
the independent variable at which dependent variable Q is beihg 
considered, and m denotes different dependent variables. Corrcsnon- 
dence between m and different dependent variable is done in a 
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systeiii8ttic way. With the order in which equations are written 

in ^1- 2) program assigns 

G = Q(n, 1) 

G' = Q(n, 2) 

G" = Q(n, 3) 

G"' = Q(n, 4) 
and T = Q(n, 5) 

T' = Q(n, 6) 

T" = Q(n, 7) 

Such correspondence in general program, is obtained by 
introducing some artificial integer variables whose values can 
be fed in the program as data. Values of these variables are 
different for different set of equations. 

The artificial variable are JJEQ, JJCI, JJCP, NKB/EP, 

LCH, RIB, ra, MI, RRV. 

In addition to these variables there are others like 
SST, PM, RR, H, BB, JA. 

Some of these variables are double subscripted, right hand 
side subscript showing the set of equations to i^'hicb these varia- 
bles belong. 

Row a short description of these variables are given b^low,, 

JJEQ is a double subscripted variable which is given 
value 0 or 1 . I f 

■ JJEQ (n, m) = 1 

it denotes highest derivative of dependent variable o.ccurring 
in one of the equations of m'^'’^ set i's the n"^^ dependent variable 
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Hence corresponding to equations (42 ) 

JJEQ (4,2) = 1, JJIQ (7, 2) = 1 
JJI.1Q (4, 2) denotes occurance of G"' and JJEOX?, 2) denotes that 
of T ” in the second set of equations. 

JJCI is a double subscripted ''lariable which tolls the 
program what ■ the missing starting values arc, !■ ike JlfQ, JJCX 
is given values 0 or 1. Now if 

JJCI (n, m) =? 1 

it denotes initial value of the n"^^ variable is missing' in 

set of equations. Hence, for equation (42), we have 

JJCI (3,2) = 1 

JJCI (6,2) = 1 

These correspond to G" and T' whose initial values are not Icno':: 

JJCE is also a double subscripted variable denoting the 
variables which have prescribed values at the end point. If 
JJCF (n, m) = 1 it denotes n'^'*^ variable of the set of equa- 
tion, has a prescribed value at the end point. For eq. {,4 

we have 

JJCF (2,2) =1 
JJCF (5,2) = 1 

These variables take: the value 1 under the conc’-'tions mentionc 
above, otherwise they take the value zero. 

LiaH is a subscripted variables. .It can take the value 
0, 1, 2, 3. If ICH (m) equal to zero, set of equation will 
be calculated as will be explained while describing subroutine B 
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If ICH (m) equal to 1 , the set will he calculated by PDM, 

If LCH(in) equal to 2, the set will be solved directly after 

obtaining approximate starting values b3?' PDM, If LCD (m) equal 
to 3, set will be solved directly and values obtained in sol* 

ing the previous set of equations will also be used in solving 
this set. 

NIB (m) and NIB (m) gives respectively the number of 
initial and final boundary condition given for the m.'^^'^ set of 
equation. Here, their values are 3 and 2 respectively for abl 
sets of equations. 

NKMP is the total number of sets of equations to be 
solved. Hence here NKMP is equal to 5 (s-s three sets of equatio: 
are to be solved here). 

With these artificial variables, program, reads the identf 
of the equation. Other variables like H,NFV,NR,SST ,EEN supply 
the actual data like range of independent variable^ values given 
at the boundary points, etc. 

A short description of these variables are given below 

,NN?(m) denotes the number of dependent variable as will 
be interpretated by the program in the set of equation, i.e. 
it is equal to the sum of the maximum order of derivatives of 

different dependent variable occurring in the m'^^ set of equa.- 

"til 

tions and number of actual dependent variable in the m set of 
equation. 
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Here for the second set of equations 

maximum order of the derivative of G = 3 
maximum order of the derivative of T- 2 
no. of dependent variahles{G and T) is 2. 

Hence 

MV (2) =3 + 2 + 2 = 7. 

HR = No. of steps into which domain of independent 
variable is divided. 

H = step size of the independent variable. 

SSI and RM are double subscripted variables corresponding to the 
boundsry conditions given at initial and end points. Here GCo), 
G'(0)5 T(0) a3?e given for the second set of equations we assigns 
these values to SST(1, 2) SST (25 2), SST(3, 2) respectively. 
Similarly the end point G" (°o) and T (°°) are given, hence we 
assign these values to FRHCI, 2) ahd FIN (2, 2) respectively. 

First initial value corresponding to the order in which 
different dependent ( variables are arranged (as interpreted by 
program) is assigned to SSI C^- , 2 ), the second to SST (2, 2) 
and so on. The order in which the dependent variable, is arranged 

G, G’, G", G'”. T, T’ , T” 

Now first initial condition corresponds to G, second to G' , third 

SST(1, 2) = G(0) =0 
SST(2, 2) = G' (,0)= 0 
S3T(3, 2) = T' (0)= 0, 


to T. Hence 
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Similarly, 

M(1) = G-" (~) = 0 
™(2) = T («) =0 

Vrrir' le IJ controls the accuracy to which boundary conditions to 
be satisfied. If iMI ia given a vali'-- 3, min: mum upto third place 
of decimal, boundary condition will be satisfied. 

5 . 3 Subroutines 

There are eight subroutines in this program. They are as 
given below. 

MIST -This subroutine along with subroutines RKM solves the 
equation by PDM. This is effectet .Jor the second set of equatio 
by putting LCH(2) =1 

AMT - This subroutine using the values of JJEQ, JJCI, JJGE, etc. 
defines some other variables like JQ, JO, letc. Values assigned 
to these variables are used in the main program and subroutine 
mCM. 

HIX - This subroutine along with BKM, solves the original equa- 
tion directly. Here original equation is the third set. Hence, 
we equate ICH (3) to two. IXie to this value assigned to ICH 
for the third set this set will be solved directly utilizing 
starting initial condition as obtained solving second set by PPM, 
RKM r- This is the most important subroutine which solves a parti- 
cular set of equatiqns depending on the value attained by KT/iP in 
the main program, IIP denotes the number o-f' the- set ? to be solved. 
When in the main program KMP attains the value 2', 
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second set of equation is solved dy the program. As mentioned 
earlier, this subroutine iterates the result to the desired 
accuracy controlled by the value assigned to M . 

AlEQ - This is a subroutine which soD.ves linear algebraic equatior 
This subroutine is called from EKM for gettirg new interpolated 
initial values,. 

TICK - This is a subroutine which calculate derivatives one order 
higher. In EKM, this is caKled to obtain higher derivative from 
lower ones. 

EUNSON - This is the subroutine where different sets of equation 
are fed in the following way. let us consider second set of equa- 
tions which can be written as 

G"' = - ( ff’ + (s-f'^)p+B(gf'H-fG''+(T-2f'GOp)-f"+G"(l-B)) 

T" = - (fs' Pr +B G s' Pr + BfT’ Pr + T ' - BT'-s') 

Now this can be written in terms of subscripted variables as 

EOMJLA = - (p(1)*P(2)+(P(5)-P(2)*F(2))*p+B*(Q(l)*P(3)+P(l)*e 
+ (Q(5) ^ 2*F(2)*Q(2)) * p - F(3)+Q(5)(1“B) ) 
where G"' is replaced by a variable FOMIlA 

variable F corresponds to f and s and Q corresponds to G and x, 
subscript of F and Q are decided by the same rule as mentioned in 
case of subscript m of double subsor-’ pted variable F(n, m). Yvliile 
writing down equations in FUNSON subroutine we should remember 
variable like G, T, G' etc. which are actually being evaluated 
in EKM by Runge-Kutta Method is denoted by Q and other variables 
like f, f ' , s etc. which has already been calculated and is being 
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utilised for calculating G-, T, G’ etc, or Q are denoted "by E* 

The statement hearing foimula given above is assigned a statement 
number given by 

I = KMP*100+M 

where M is eq.ual to the position of a equation in the KMP'^^ set 
of equations. The above equation is the first equation of second 
set. Hence, 

N = 2*100+1 = 201 

Hence, equation (44) is transferred to EUHSOH subroutine as 
( for p = 1 ) . 

2ri lONULA = _ (P(1 )*e(2)+E(5)-E(2)*E(2))+B*(G(1)*E(3)+E(1)*Q(3) 

+ (Q(5)-2*P(2)*Q(2))*1-F(3)+Q(3)*(1-B)) 

Similarly, equation (45) is transferred as (for Pr = 1 ) 

202 JOmhk = - (E(1)*F(6)+B*Q(1)*P(6)+B+P(1)*Q(6)+Q(6) 

-B*Q(6) - P(6)) 

When the third set 'is to be- solved, corresnonding equations 
are transferred ' to the FUNSOIT . subroutine as (for p = 1, 

Pr = 1) 

301 ..FOMUL/K = - (Q(1 )*Q(3)+Q(2)**2)) 

502 FOWUL/v = - Q(1)*Q(6) 

BST - If solution for the initial value of the parameter can be 
obtained directly as in the ease of equation (41) where solution 
is obtainable in a closed form as given by (41a), then the formulc 
giving the solution for the initihl value of parameter is fed in 
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the subroutine BST. This is ochioved by giving the value equal to 
zero to the variable ICH for the set of the equation for which the 
solution is obtainable in a closed form. Here equation (41) is the 
first set. Hence v/e put ' 

LCH (1) = 0. 

5.4 Some comments on the General Program? 

To use the general program, most economically and effective!' 
it is suggested that first the problem should be solved by PPM 
taking large 'step size of the perametur. After obtaining the 
solution at the desired value of the parameter, we can solve the 
original non-linear equations directly with starting values obtainei 
oy PPM in the first step. As the assumed initial values are very 
near the exact valu"’, solution will converge to the true solution 
of the desired accuracy in a few 'x-terations. Solution obtained by 
PPM by taking a step size of p equal to .25 and the solution obtain* 
by solving the parent equations directly utilizing the initial boun* 
ary valu® as obtained by PPM, are given in a tabular form on the 
next page. 

In the main prcrram, double subscript variable AF always 
denotes the variable whoee values are to evaluated like (f, f ' , s,s 
etc. The other double subscripted variable is Q which is the vari- 
able whose values are known. In case of PPM, Q denotes variables 
like (df/d:p), (df/dB) etc. 
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Table 1(a) 


ETk 

f 

s 

O 

« 

o 

0.000 

1 . 5000 

0. 5 

0.1552 

1.352 

1.0 

0.5171 

1.217 

1.5 

0.9722 

1.113 

« 

O 

1.460 

1.048 

2.5 

1.936 

1.017 

3 . 0 

2.453 

1.003 

3.5 

2.950 

1.001 

4.0 

3.448 

1.001 

4.5 

3.946 

1.001 

5.0 

4-445 

1.001 

5.5 

4.944 

1.000 

6.0 

5.443 

1.000 

6.5 

5.942 

1.000 

7.0 

6 . 442 

1.000 

7.5 

6.942 

1.000 


T Solution obtained, by solving 
equation (42) by PDM taking 
B = 0.25 for p=0,?r=1 


Table 1(b) 


ETA(y-|) 

f 

s 

0.0 

0.000 

1.500 

0.5 

0. 1565 

1.352 

1.0 

0.5213 

1.217 

1.5 

0.9800 

. 1.112 

2.0 

1.472 

1.048 

2.5 

1.971 

1.017 

3.0 

2,472 

1.005 

3.5 

2.972 

1.001 

4.0 

3.472 

1.000 

4.5 

3.972 

1.000 

5.0 

4.472 

1.000 

5.5 

4.972 

1 , 000 

6 . 0 

5 . 472 

1 . 000 

6.5 ■ 

5.972 

1.000 

7.0 

6 . 472 

1.000 

7.5 

6.972 

1 . 000 


Solution obtained by directly 
solving the parent nonlinear 
equations starting value of whic> 
is obtained by solving equation 
(42) by PDM. Here also p=0,Pr=1 . 
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The program can also he utilised to solve more than one 
set of equations where results of one set of equations is to he 
used in tKei solution of the next set as is required in solving 
second order boundary layer equations where values from the first 
order is required to solve the second order equations. In this 
case Q can he given values corresponding to the first order 
solutions. This is achieved in the program hy giving the value 
3 to the variable LCH for the set of equation corresponding to 
second order boundary equations. 

General program is given in the Appendix 3. This program 
can solve a set of equations containing 25 (computer) dependent 
variables. On IBII 7044 computer, some 4500 memory cores remain 
unused. Hence by changing dimension, this program, can handle 
equations with 35 dependent ■''"ariable. Maximum number of boundary 
condition that can be satisfied at the end point is 10. Time take 
by this program on IBM 7044 to solve equation (42) by PBM and then 
to solve equation (43) directly is 6 minutes 53 seconds. 



DISCUSSION 


Results obtained for Pra.ndtl nurr'ber equal to unity tally 
quite fairly with We, given by Cohen and Reshotbo (Ref. 4). 

Gene rv;l. trend, of the reeult ip the same e.s siven in fVt>^t reference 
(4). haxiiTium deviation of our result frori that given in that 
reference is ,5 percents This difference ca-n still be reduced 
by takina shorter step size of paiaioeter p and thaat of the 
independent variable . Prom the result obtained .for Prandtl 
number eual to unitv, we can hope results obtained for other 
Prandtl number are also- reliable so far a.s the solution of 
equation (15) and (16) is concerned. 

For Prandtl number not equal to unity, solution o'’'>t3ined. 
in general is the complementarv function of equation (13) hut 
gives the near exact solution for low mach number flows. For 
high mach number (M - <») right hand side of eqnation (13) 
approaches a constant value equal to two. , With slight modifi- 
cation in programme results can he obtained for high mach number. 
Thus solution obtained for Prandtl number not equal to one may 
not give exact quantitative value in all cases. Nonetheless, 
values were calculated for Prandtl number not equal to unit^!" 
as it gives a qualitative understanding of boundarv layer. 

Our main aim was to solve the compressible bo-undari,^ 
layer equation by parametric different is-t ion method and we show 
that this method can he effectively auplied to solve the prohleru 
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APPENDIX -i 



l-MLB I 

liFl'A = 0.0000 PRANDTLE NO. = 1.0000 

V..:!.(:nirY GRAOIKHT at wall = 0.4696 

Tf-MPFRATURO GRADIENT AT WALL = 0.4695 


r I SPL ACr.PhNT THICKNESS = 

-0. 

0001 


MOMEMTtjM 

TUICKNE 

SS ■= 0 

,4696 



THE RMAL 

TMICKNES 

s = ' -1. 

2169 



REYNOLD 

ANALOGY 

PARAMETER = 

2.0004 


ETA 


F 


DF 


0. 


0, 


0.l428e- 

07 

' w J ' ' 

j 

0,5064!=- 

CL 

G .2342 F 

00 

-■ . ] -r: ; , 

i 

0.2330c 

00 

0.4606E 

00 

, l o: 

1 

0.5150E 

0 0 

0.66i5S 

OO 


C! 

O.oaOciE 

00 

0 .6 16 7 h 

00 


Ci 

0. 13225 

01 

0,91 68 F 

00 

. .Kor’'' 

: Li 

0.1796!:: 

01 

0.9 69 1C 

00 


■' Ci 

0. 2286E 

01 

0.9907C 

00 

' :: . A (': ( ) • :■ ' 

f t 

0.27P4E 

01 

0.9978E 

00 

‘fj . A ‘"i A ' 

' DI 

C.32P3E 

0 1 

0.9996E 

00 

0 . t ' 

c; 1 

0.37B3E 

01 

0.9999E 

00 

J. OOOtA 

i. 

0.42B3E 

01 

O.IOOOF 

01 

0. (/COD': 

: (;1 

0.47P,3E 

Oi 

O.IOOOE 

01 

0 .£ 'VD-I ' 

'' i.l 

0.5283E 

01 

O.IOOOE 

01 


i: Cl 

0.5783E 

01 

O.IOOOE 

01 

n. ?5c fv 

t Cl 

G.62B3E 

01 

O.IOOOE 

01 


SH 


0 . 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


o.oooi; 


s 

. 2 342 E 00 
.4606E 00 
.6614F 00 
,8167E 00 
.9i6&E 00 
.96 OOF: 00 
.9907E 00 
.9978E 00 
.9996E 00 
.9999E 00 
i.lOOOE 01 
f. IDOOE G1 
I.lOOOE 01 
I.lOOOE Cl 
I.lOOOE 01 


ta&le X. 


BETA = 1,0000 PRANDTLE NO. = 1.0000 SW = 0.0000 

Vri.nriry gbaoie-.nt at wall = 0.648,6 

Tl:: ' riJl-T GRAOIHF^T AT WALL = 0.5065 

I) i StM. Ar.r«E-jT THICKNESS =■ -0.1575^ 

MOBFNTUB THICKNESS = 0.4036 

T!if'B"AL THICKNESS = -1.1374 


evr^iCLf) 

A IV A LOGY 

PARAMETER = 

2.5612 




ETA 


F 


OF 


s 


f; * 


^0:. 


B.l « 12 ■ 1' '.1 

07 

0. 


, *'Br ;/'■■■ 

'i 

0.79H7t- 

01 

G. 3149E 

CO 

0.2524e 

00 


! 

C . 3 0 6 7 E 

CO 

0.5B17E 

CO 

0.4936E 

00 


i. i 

G.6A9'^E 

00 

0. 7762E 

00 

0.6998E 

00 

'■ . / 0. '■ . 1 

01 

C;.1071E 

Cl 

0.8 97 IE 

00 

0.8492E 

00 

1 

0] 

0.151TE 

0 1 

C.96(}?e 

00 

0.9374E 

00 

f r ‘I \ C' 

' . •# . 

0 1 

C,2U25E 

01 

0.9873E 

00 

0.9789E 

00 


Cl 

0.2521E 

01 

0,9967E 

00 

0.9943E 

00 

^ A ' 

1 

0. 2020E 

01 

0.9993E 

CO 

0.998£E 

DO 

0 * A 

Ci 

0 . 3 5 2 C £ 

01 

0.9999E 

CO 

0.999aE 

00 

i' IB; ! MB f : 

0,1 

0.4020E 

01 

0. lOOOE 

01 

0. LCOOE 

01 

f ; * E’) (, ?' 

(. l 

0,45?OE 

01 

0. lOOOE 

01 

0. lOOOE 

01 

H .ocoon 

C1 

0.5020E 

01 

0. lOOOE 

01 

O.iOOOE 

01 

6 0-,;C:!,, 

1 1, 

0.55?oe 

01 

0 . loooe 

01 

C.IOOOE 

01 

B.TCOBE 

■ Cl 

0 * 6020E 

01 

0. lOOOE 

01 

O.IOOOE 

01 

0. TOOOu 

. -Cl 

0.6520E 

01 

C.IOOOE 

01 

o.ioooe 

01 
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TA&te 3 

RFTA = '■=.U0-')0 PRANDTLE MO. = 1.0000 SW = 0.0000 

■/'IfOlTY 0)lAnir-NT AT WALL = 0.7384 

Th ' ’I'i' \ I I.W-ir GRADIENT AT WALL = .0.5205 


s^m.acf: 

i'0-:.jr THICKNESS = 

-0. 

2060 




M c f.,; j 

IHICK.HE 

SS =, 0 

.3839 





r RMAL 

THICKNESS = -1. 

1109 





YMCI f) 

analogy 

PARAMETER = 

2.8374 




'E T A 


F 


OF 


s 




0. 


-0.3332E- 

•07 

0. 


r“ /*'!, f\ 

CO 

O.OOOIE- 


0.351TE 

00 

0.2593E 

00 

T 

.«» i " 

i'; 1 

0. 3385E 

7*^ f , 

0.6270E 

00 

0. 5056 E 

on 

1. 1 ''V; F". 

l 

0 . 70 1 4 E 

0 0 

0. 3n99E 

00 

G.7131E 

00 

^ “? ' F ' 

*' 1 

G. 1 I36E 

0.1 

0.0156E 

00 

0.8596E 

00 


(} ] 

0. ] 608 E 

01 

0.9682E 

GO 

0.9434E 

00 

. :■ i'’. ■’ ' ‘ ' •' 

TFI. 

0.2099E 

01 

0.9900E 

00 

0.9F15E 

00 


Cl 

0,2596e 

0 L 

C.9975E 

00 

0.9952E 

00 

. 41 -im'; 

Cl ' 

C.3095E 

01 

0.99956 

00 

0.9990E 

00 

, !'4 1 

F 1 

0.3593E 

01 

0.9999E 

00 

0.9998E 

00 

♦ ’"/ Vj, 1,' "' 

Cl, 

0.4095S . 

01 

0. IQGOE 

01 

0 • lOOOE 

01 

'■'1 '‘'"j / ■' ! 1 

C 1 

0 . 4 5 9 5 E 

01 

0. lOOOC 

01 

0. lOOOE 

01 

. 6000 = 

u 1. 

0.5095E 

01 

O.IOOOE 

01 

O.IOOOE 

01 

. OOOO'i,' 

Cl 

0.559‘5E 

01 

O.IOOOE 

01 

0. lOOOE 

01 

. /oo='r 

r: 1 

0.6095E 

01 

O.IOOOE 

01 

O.IOOOE 

01 

. 1 i i C f'-,' 

(M 

0 . 6 5 9 5 E 

01 

O.IOOOE 

01 

0. lOOOE 

01 



TABi-E 4 

['.■^TA = J.aooo PRANDTLE NO.. = l.OOOu SW = 0.8000 

Vn.00lTY LrlAniFMT AT WALL = 0.A696 

Ti::-^P'--PATljHS-- GRADIENT AT WALL = 0.0939 


DI sRLAr,rf’'"Mr thickness = 0.9734 

MOMFMTli!'": THICKNESS = 0.4696 


rm-R'MAi. 

THICKNES 

S = -0. 

2434 





Ri:yNCL.D 

ANALOGY 

PARAMETER = 

2.0004 




&TA 


F 


OF 


s 


0 . 


0. 


0.14286- 

07 

0.8000E 

00 

. S c ■.-I ;• r 

C 0 

G.5864E- 

•01 

0.2342E 

DO , 

0.84686 

GO 

V. I oo:,:'- 

■ 1 

0.2330E 

0 0 

0 .4606 F 

uO 

0.89216 

00 

v: . 1 

1 

0 . 5 1 3 0 E 

00 

0.6615E 

00 

0.93236 

00 

0. 

1 

0.8868E 

no 

. D.e.l67E 

00 

0.9633E 

00 


.. l 

0. 13?2E 

01 

0.9L68E 

00 

0.98346 

00 

’ ' % T ' ■ ..J . ' 

'' 1 

0. 1,796E 

Oi 

0.96915 

00 

0.99386 

00 

, '0'':.>V,r; 

Cl 

0. 2286E 

0 1 

0.9907E • 

00 

0.99816 

00 

AC .rOE 

■i 1 

0.27B4E 

01 

D.997BF 

00 

0.9996E 

00 

; . . A CC'Ci' 

1 

Q.3283E 

01, 

0.9996E 

00 

0.99996 

00 

. : . 3 C L- A’ 1 ; 

•:.,i 

0.37B3E 

Oi 

0.9999E 

00 

0. 10006 

Cl 

' • 4 ‘i -j CT'., 

01 

0.4283E 

01 

0. lOOOE 

01 

O.IOOOE 

01 

u. 6 OR Of" 

Cl 

0.4783E 

01 

0. lOOOE 

01 

O.IOOOE 

01 

o.ooQCi:: 

Cl 

0.52836 

01 

O.IOOOE 

01 

O.IOOOE 

01 

TOO'-'t: 

Cl 

0.5T83E 

01 

0. lOOOE 

01 

0.10006 

01 

.■..TSoiO, 

01 

0.6283E 

01 

O.IOOOE 

01 

O.IOOOE 

01 



RfTA = 

TAgLe 

1.0000 PRANDTLE NO 

s~ 

. = 1.0000 

SW = 0. 

3S 

,8000 

Virl.nc I TY 

ORADIfjNT 'AT WALL = 

1.1234 




TCFPFR AT 

UHE GRADIENT AT 

WALL = 

0.1119 




DISPLACE 

HENT THICKNESS = 

^ 0.4950 




M OHENTUM 

THICKNESS = C 

'.3163 





THERMAL 

THICKNESS =. -0. 

2088 





REYNOLD 

ANALOGY PARAMETER = 

4.0165 




fc T A 

F 


OF 


s 


0. 

0. 


-0.9810E- 

07 

0.8000E 

00 

0.5'CO'Oe 

CO 0.123iiE 

cc 

0 *4622E 

00 

0.a556E 

00 

D. ICE or: 

01 0.4320E 

GO ■ 

O.T442fc 

00 

0.907&E 

00 

w. 1 

il G.845BE 

on 

0.a934E 

00 

0.9494E 

00 

? C T‘ VC 

Cl n.l312E 

01 ' 

0.9616E 

00 

0.9770E 

00 

A • R f 

Cl 0.1801E 

01 

0.9881E 

00 

0.9914E 

00 


Cl 0.22'17E 

01 

0.9969E 

00 

0.9974E 

00 

'3 

Cl 0.2796E 

01 

0.9993E 

00 

0.9994E 

00 

Ir # A t '. 0 f'i !'■" 

Cl rj.3296F 

01. 

0.9999H 

00 

0.99996 

00 

A '3e> v,iF 

Cl 0.3796E 

01 

O.iOOOE 

01 

O.IOOOE 

01 

I'i # o> C u f' ' 
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GFTA = ?.00n0 PfUNDTLE NO. = 1.0000 . SW = 0.8000 

V'-l.r-rifY ORAtM-NT AT WALL = l'.5127 
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BKTA = 1.0000 PR&NOTLE NO. = 1.0000 SW ■= 

VMJJCITY GOAfMENT AT WALL = 1.A904 
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"•i ''V 
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O.IOOOE 01 
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O.IOOOE 01 

O.IOOOE 01 

' j « '*,1 0 0 

r ,01 -.O.SQAAE 01. 

O.IOOOE 01 

O.IOOOE ^01 

0 . 6n:]G 

0- 01 01 

O.IOOOE 01 

O.IOOOE 01 

G . 6 A ) 

f' 01 0.6044K 01 

O.IOOOE 01 

O.IOOOE 01 

0 . 7CG0 

oL ' 0.6544E -01 . 

O.IOOOE 01 

O.IOOOE 01 

0 . ?Gv'0 

.'31 0,,?044e 01 

O.IOOOE 01 

O.IOOOE 01 



r-riA ^ PR&WOTLE.MO. = 1,0000 SW = 0.8000 

Vi 5 if'! I'/ ' 1 r NT AT WALL = 0.3952' 

It it’i M;;.- OLAOTTmt at wall = 0.0908 

1 1 r, tM 0 1 i T ri' ■ A ^ c ■. * 1 . 071 8 


'Hf I ii nur.KMFSS = 0.4929 

'Tl-' T!'!rK?!!::i,S = -0'.2'50? 
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1 
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f i . 1 2 c 
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: ■■■■i 

!>, I695t 

(j i 
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r,.p]J>2’C 

01 


f 

0.2 6 7 91-3 

(,) 1 


i, L- 1. 

0..3178F 

01 

: • ■■ ' ’ . ’ ; • ‘ 

’ i 

0, 3673E 

01 

i" f. , - 

(' A 1 

0.41 7 7E 

01 

/ 1 ■ ■ ' S ' 

i. 1 

0.46776 

01 

« 


0.'V17 7F 

01 


1 

0.567 76 

Cl 

''PSU"i 

t; 

0*6 1770 

01 


1.T413 

OF ■ S 


0.6726E-09 

0 » o'O D vJ E 

GO 

C.2C35E -00 

0 • 8 A’ 5 3 E 

0 ' J 

G.4132t 00 

0.8S93E 

00 

C.6123F 00 

0.9289F 

DO 

0.77766 00 

' 0 « 9 6 C 3 E 

00 

0.R926E 00 

0.9B13E 

GO 

0.95 746 OC 

0.9927E 

00 

0.9B63E 00 

D.9977E 

00 

0.9965F 00 

0.9994E 

00 

0.9993E GO 

0,9999E 

o'o 

0.9999fc 00 

O.IOOOE 

01 

O.lOOOE 01 

O.IOOOE 

01 

O.IOOOE 01 

O.IOOOE 

01 

O.IOOOE 01 

O.IOOOE 

01 

O.IOOOE 01 

0. lOOOE 

01 

O.ioOQE 01 

o.inooE 

01 



TAei-/' ^7 

is. >■ T A - ~ . I ? f) n A M !') TIC ;M'0 *=,1.0000 SW= 0.8000 

V i .; !i'-' f/{i0>is-0r AT WALL = 0.3049 

Tf:Hi:>rr’Al i„i;0,-: ORAOIFNT . AT WALL = 0.0865 

r'- 1 i’l. A {,: ■ m f.i t i i f r. ' a ^-t. s = i . 2 1 4 1 

I h: . f.noK^lFSS = 0.5217 

Tf.iCKNrSS = : -0.,.2617 


Yf4t::l 1' 

1 A, A LOGY 

li /■-, l\ /'^ !/: A r ^ f; 3= 

r 

1.4102 

■Of 

S 



■*" 1 .) ♦ 

'i « 4 ■' 1 1 5 F “ 0 1 

-0M019E-07 

G.164GE 00 

O.faAO^AH wO 

0.B432E 00 


, ■ “t 

A. 16845 DO 

0.3506E 00 

0 . 3 fc 5 3 E 0 C) 

«> i 

' ' ^ 1 

a.oopoe go 

0 . 5 4 3 1 E G 

0.V240E 00 

• j ; • 


C. 70 86 8 OG 

C.7103F 00 

0.9557E 00 

■ <a / " ' ' 

;' . . ^ 

o.rioiF Cl 

0.35260 0 0 

0.9780E 00 



0.15538 Cl 

U.0 363E OC) 

0.9909E 00 



.■1.203 IE 01 

0.9T77E 00 

0.9970E 00 

-.Ai . 

' 'I 

U.2526fc 01 

0.9938E 00 

0.99928: 00 

^ A^. : 

: A J- :A ] 

0.3025H 01 

0.9986E 00 

0.9998E 00 


M'iA Cl 

. 0.3524e 0,1 

0 .999 B L 0 0 

O.IOOOE 01 


^ 1 ' r 1/ i' . . ' 

. a, 4 024 E/,01 

O.IOOOE 01 

O.IOOOE 01 

» C 0 ' 

,f ’.A 

0.4 5241: 01 

O.IOOOE 01 

O.IOOOE 01 

u * ' 

,4.7 ui \ 

C.5024E Oi 

O.IOOOE 01 

O.IOOOE 01 

* 7 t"! 

40 it 01 

,0.5524E CM , 

o.icooe- 01 

O.IOOOE 01 


4'V''" ■ CCi 

A 0,.6024e CM 

0. LQOOE 01 

O.IOOOE 01 




TA filE 



; > 1' A - 

- » . i ny’i pMANDTLr nh 

. = 1.0000 

sw ■= 0.8 000 

VM I 1 '' 

GR AO I ENT' AT WALL = 

0,1794 


; A ’ 

SiOr GRADIENT AT ^WALL = 

0.0792 


OISf'l'^CL 

»'.;MT TH if.KGCRS ^ l./tTlT 


f’j ,,v;' !■: III?' 

liilCKDFSS = 0.!3599 



! )i*. a;-'-' 1 

TRICK ME VS -0.2:.?, 6 



rfymi,..! r- 

AKAIOCY PA,'. A NET PK - 

0.9055 


■ r f. 

" ■ F "" 

OF 

s 


"0. 

-0.1273E-07 

0.8000E CO 

* ■ ' ' 

. , 0.AG57E-G1 

C.1086E 00 

0.8396E 00 

«« i 

' ! ':>.11A7T OG 

0.2537F OC 

0.8TS6E 00 

■ . 1 - ■■ 

' 1 O.RHACE 00 

0.4267E 00 

0.9153E 00 


'■! O.GApTF 00 

0.6075t 00 

0.94706 00 


. 1 0.8 33 IE CO 

G.76B1E 00 

0.9712E OC 

■s ■' ‘ 

^ 1. 0.130^!r 01 

0.8854E 00 

0.98686 00 


C. 176 OF 01 

0.9 538E !jO 

0.9S5GE CO 

■ ^ ^ 

0.2251,0 01 

0.9850E 00 

0.99856 00 

■**'*>' V, . ■ ■ ■) 

■■■■1 0.?74,TE ot 

0.99616 00 

0.9996E 00 

^ A,: 

'■■ . n*.32A6t' 01 

0.9992E 00 

0.99996 00 

= .? Hf *V '.H' ' ■' 

; i t, 0-.374,SE 01 

■ 0.99996 00 

0.10006 01 


. •■,| Q»4247),E 01 

O.IOOOE 01 

O.IOOOE 01 

. ^ i'. f 

r , , 0.4r46e 01 

0.1 CODE Ul 

0.10006 01 

A' . 7 i . 0 C 

r ;,i O', S'? 4 66- 01 

O.IOOOE 01 

O.IOOOE 01 

A, ?AOO 

oi, :G4,5746E O'l', 

O.IOOOE 01 

O.IOOOE 01 



' 2 .^ 

= 0.3000 SW = 0.8G00 

0.4149 
= 0.0597 


0.9015 


■TAgte 

BETA = -0.0500 PRANOTLE NO 

VELCCfTY ORAPIF-NT AT WALL = 
Ts-.YPORATl.Rh GRADIENT AT WALL 
niEPl.AOEPEMT Ti-:i~7 
POf-^EMTUM THICKNESS = 0.4840 

THERyAI. THICKNESS = -0.3866 

REYNOLD ANALOGY PARAMETER = 


ETA 


F 


0- 


-u . 


■j. LCOO " 
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•01 
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Ci 

C.2i26F 

Ou 
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Cl 

0.4774E 

00 
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G 1 

C . 8 3 3 9 F 

GO 
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ri 
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0 1 
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Cl 

G.1727E 

01 

0. 3 5(jO- 

01 

0.2216E 

0 1 

0.4 Cl: or 

C I 

0.271 3F 

0 1 

0.4 5001- 

Cl 

0.3212E 

01 

0. OCOOr 

Cl , 

0.3712E 

01 

0. 5‘.^COE 

Cl 

0 . 4 2 1 2 F 

01. 

0.6■OOOf: 

Cl 

0.47 12 r 

01 
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C 1 
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01 
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Cl 
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01 
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Cl 
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OF S 
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,0 , S : 0 5 E 

00 
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00 
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CO 
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GO 

.0.85940 

00 
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00 
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00 

0.7908E 

00 

0.9146E 

00 
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00 

0. 93800 

00 

0.9622E 

00 

0:9575E 

00 

0.9865E 

no 
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00 

0.9973E 

00 
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00 
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00 
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01 

O.IOOOE 
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"TA6ljf 3X> 

BETA = -0.0625 PRANDTLE NO. = 1.0000 SW = l.OOOO 


ViJlB'.ITY GRAiOIBMT AT WALL = 0.3B14 

TE-'- 'tA H:Rr OBMJTtNT AT WALL = -0.0000 

fl! ','^1 ACO THICKNESS = 1.3407 

HOMEMTUN THICKHESS = 0.4962 

THEKf^AL THICKNESS = -0.0 COO 
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SETA = -0.0625 PRANOTLE: NO-. = 9.0000 SW = 0.8000 

\/i'LCC.lTY 'Vl.MMf'MT AT WALL = 0.3888 

T £ P E R A T IJ R E 0 R A D I 'lAJ T A T W A L L ' = 0,1889 
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thickness = 0.5279 
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f-EG'-'''' FrKTR-'/; SOURCE LIST 

’ 5^4 sni, 'KCr ST■VTE^*F^iT ■ 

0 t-inric . ■ ' ' , , ■ 

1 I ] )• cNSIij I f ( ?10) ,Fl(2in) ,F? (210) ,F3 (?10) , S { 210) , S L { 210) ,S2 ( 21 

1 -OU ::>1C) ,G1 { ?10) ,G2 ( ?10) ,G3(210 ) ,T1 ( 210) ,T2( 210) ,T3( 210) , 

2 AhTA{?10) ,DF{210) ,BS(2ir;) 

3 ,r '(o! ,i):. (6) t AA(b) ,i:iF{2in) ,AFP(?16) ,SWW(6) 

2 AL lao, 

3 4 ri'Ki'\'‘.T(/.'?:;X,2 7!4y'=LCCr"Y GRADIENT AT WALL =,F9.4) 

2 F::jRi'AT(/2'X, SOHTFYA-'rXATURE GRAOIEHT AT WALL =,F9.4) 

3 iCi -i;FFi\T(/?5X,3HeTA,12XtlHF,14X,2!^DF? t3X,lHS/2X) 

6 F ' T( /2uX,24hDI 5PLAFFMF"iT THICKNESS =,F9.4) 

7 9 FOi-!' AT( IMl ,2GXf 6HFFT3 =,^9.4 t3Xt 14H'‘RAMDTLS NO. =,F9.4 ,3X, 

] ~ t F » 4 ) 

]'■ 7 Ff;Ri-AT(/2CXr?0HN0Yc3T!;M THICKNESS =,F9.4) 

11 3 f !'3' 7 T(/?cx,19hTHCRr'Ai. THICKNESS =,F9.4) 

17 99 F.';k- A'.T{ /20X.27HRFVNCLL Ar’ALOGY PARAMETFR =,F9.4) 

I Y 11 F f :i ?./: A T ( / 2 X , F 1 1 . 4 , 4 X , H 1 1 . 4 , 4 X t F 1 1 . 4 , 4 X , e 1 1 . 4 ) 

)4 P.FAilGFO ,H,DP 

I H '414 F(1'H’AT{4E15.5) 

If'. 999 FuR:<AT(F4.2,F4.2) • 

17 ??7 FnKpAT(QEl?.4 ) 

20 55 9 FriYAKFS.?, lOX, r IP . 5 1 lOX , E ] 9 . 5 , lOX , I 5 ) 

21 41,9 f-ORf-‘AT(3215.5) 

??. CALL FLUN(29C0) 

23 KST=l 

24 SV.'H{1) = 0.2 

2‘^- SHWi2)'="C,6 

26 ' SyWS3)=l.6 

27 SWw(4)=i. 

30 SwW(5)=?. 

31 IKK =4 ■ 

32 PR=5. 

3:L 147 n0331YST=3, IKK ' ' . 

34 JA=21 

33 ' .9A.-191 

36. , ■ A'=l. 

37- ' ■ ASTA( 1).= 0. 

- 40 . .; ,001.13,1 = 1 ,2 ^ . : 

,/..7,:4t/ \.A''9... :,:'S-=o.' x .. x - ' . , , ' , . ■ , , x i, 

.,y7.:'^:3 019,0.' : . ■ • . ' ' , . . ' ,. ■ ' 

, :'.44 .. .. ^ . W” A,';.. .. .' 

45 F.TA=C. .. : ,■ ■ ■ ' X , ' 

' ■ 46 ' '0=1 '. . ' 

4 ■? F ( N ) ='G . ■ 

90 222 eK0'=HY'7 . : . ' ' x . , 

51 OL0=H^W , 

57 nxO=-H*(F('''i)=5‘W) .' ■■ .' 

53 0K.1=H*(Z + 0L.0/2.) ", ■ 

54 0 L 1 '=F* T ( ij?"*'!) / 2 » ) ■ ' . ,.; ' . 

55 r;x./:i=-Vi4{ {'F (N)+'DK0/2. )^'(W+0M0/:2. ) 1 7,0 7 ■ , , 

56 C K 2 .= H ’i' ( Z + iJ L 1 / 2 . ) ,. X ’ - '' ... ■ ' ■ ' X ■ 

57 OL? = H* ( W+0?41'/?. ) ' , . . .' ...'^x'-'O, 

60 f.:Y.2 = -i''XM{r{M}+0Kl/2.)4(W + 0Ml'/2.)) ..'O;. ,. 

7,1 i x< Xsii'-M /•« 0L2 ) 

6? r'L3 = H’!''( '.vffi.'HZ ) 




for TRAM SOUhCF list 

•? r ' 

i , 


sr,! CO s^OTE'-'r-MT 

6 


r,,3 = -li4( ( F ( 1) +nK2 ) " ( ' 2 +DH 2 ) ) 

{■y, 


r { 7'4 1 ) = r { M) +1 . /6, =!.■ ( OK : +? , *DKi + 2 .4'nK.2 + 0K3 ) 

l"; ',i 


Z 7 ■•■ 1 . / (i> . * ( D L • ) + ? , 4 D L 1 2 . 4 r; L 2 + 0 L 3 ) 

() 


= 1 + 1 . / 6 , =('• ( 0 17 0 + 2 . 4 0 r" 1 + * 1 y ' 2 + n ?■•■' ■•i ) 

hi 


ir 4)113,114,114 

)’(■'. 

.i 1 -.i 

'TA=r TA+H 

1 1 


f.+-l 

/? 


C/' '^C >?2 

7 3 

1 1 4 


7 4 

1 1 '* 

r T‘, , T T MM r 

0.- t j, ’ '■« l,> 

!(■' 


CALL’’' ICK{H,OA, FfDIF ) 

7 / 


H = I 

i; ’) 


5 K = I. 

] '■■; J 


OiK = 1 

1 C /' 


sy=Shh(oST) 

1 i : i 



1 .y- 


S( I )=SC 

1 

r* ! 

' ) 71 

r:i r a 

17'; 6 


L = 1 

] 1/7 


i:N(1}=o. 

! 1 r.’ 


r.y ( 1 ) =n. 

i 1 1 

K n, 

A = 0 

‘ 12 


i 1;,/ 2 i I I ~ 4 f ‘i> 

i 1 3 


CALL i’’LlJM{ 3 00) 

i ,1.4 


ir- ( 1 .Gr,2)A=i./2. 

117 


IF n . 0 T. 4 )A=i. 

1 2 2 


FLLL=FtL)+DlF(L)*A4H 

1 2 3 


LiM{ i 1 ={P+0'J{ I-l 

1 ?4 


OK* { \ ) =-H4pR«FLLL* ( P + DK ( I-l ) 4A ) 

12? 

2 1 1 

COM TIN! IF 

i 2 7 


on 0=* { 1 . 76* ) * ( DM (21+2. ^-DH ( 3 ) +2 . 4DH { 4 ) + 0 V { 5 ) ) 

1 3C 


000= { 1. /6. )4 (0M( 2 ) + 2-4l)N (31+2. 4DN(4) + DN( 5 ) ) 

1.31 


SI L+1 l^SID+nnM 

1 32 


p=p+-Dr.N 

133 • 


,SS-S(L4-1,) „ 

1 34 


IF {L-NA)66,77,77 

13? 

fc h , ’■ 

1 = L4-I ^ 

134 


.GtlTOFS ■ . ■ , , 

1 37 

77 y 

IT {J0S{ SI NA)-1. ) -.00001) 72,92, 88 

140 

a8y,V 

KA,4lH,)*S{.?JA.) . 

141 


■ iMM::r-‘,2),.a9,9.G,9o 

142 

„ : 8-9 /, 

■'1 r=2''-''' ■ 

143 


w=w'+o'.:,50oe ' ■ ■ ■ 

1 4 4 . 


,MK.=NK+1 , . , ■ 

1 45 


001056 ■ , 

1^6 

9 r) 

IF (KK.0Q.2)Gt’T092 : 

,1 3 1 

91 

IF { AP S I AA (2 )-AA { 1) ) .LT. .01 )G0T089 

] 54 


L= { A A (2 )-L4 AA ( 1 ) +VMI. ) / ( AA{ 2 )-AA( 1 ) | 

I '5 5 


K,K,=KK + 1 

1 54 


0i:TQF6 ■ ; ■ 

i ■:> ? 

/ 2 

i';r444 K= i , N.A ' ■ ' . , : 

i 60 

4 4 

AL T A { K>1 ) =AET.A { Ks ) +H ■ ■ K. 




MEGOO'V 


'US^r.F ST-r 


FORTRAN SOURCE LIST 


I F IKS T. IT .3 ) 00102 

)T DO10aH=l,NA ' ■' . 

FUM)=F(M) 

S1(M)=S(M,) 

CCjNTIMUe 

» Tl, AZ, BZf LM ) 

i-oliJl !'^=|, ,MA 

F2 ( M ) =F 1 { M) +Gr( 

)1 S2 (M)=SUM) + Tl (M)*DB/2. 

C=B+0C/2. , , , 

CALLKUTTA{NA,!>R,F2,S2,CtHrG2,T2,AZ,HZ,LH3 



P 


:r?{r,)+G 

? {F:3*|}C/2. 


,. (o 
■ ] :■ 

F * • .5 

F3 ( R ) 

: : ' + 1 ' = >■ 

= S2(^'') + 

T2('M)=5=D&/?. 


] 1 

C 

l„ 1 ^ 

til TA(NA 

tPR,F3,s:-?r: 

,H,G:,T3,AZ,B 

1 ? 

2 i: 

1 i' 'a 
•L ' , , ■* 

= 1,NA 



1 

i ( 

=C 

3/6.* (G 

l(M3+4.«:-i2{ 

F )+G3(N) )+F(N 

'• 1 4 

1 S 1 

) ZZi'' 

■' /6.*( Tl {M)+4.*T2 { 

0) + T3I'‘'') )+S(0 

^ I (■:: 

IF 

( J-J 

A ) 1 OZ , L 

05,103 


^ i '7 

] :*"* /- p ^ 

! ■! 4- r ' 




'■ 2 ■ 3 


j+i 




' ? 1 

I r 

{{Kb 

T.FG. 1 ) 

. AND. ( J. EQ. 

11) )GCTQ? 

2 ^'V 

IF 

{ { F S 

T.FO. IJ 

.AND. ( J.£0. 

21) )G0T02 

' 2 r 

1, 

A” = ( 

I.F0.2) 

. OR . ( J . F y . 6 3 . OH . ( J . EO . 8 ) 


n 

{ { K S 

^ ♦ i: P? » 2, ) 

.AND. LAG 3 GOTO? 

■ 1 

IF 

( { F b 

T.FU..2) 

.AND. ( J.GT. 

8) )GnT0342 

'3A 

T F 

{ ( K S 

T.Fa.33 

.AND. ( J. tO. 

11 3 3 GOTO? 

‘4 1 

ir 

{ ( K S 

T.rg.S) 

•AMO. ( J.EO. 

21) 3GOTC2 


IF 

{ ( K S 

r.F</,4 ) 

.AND. ( J.EO. 

113 )G0T02 

4 1 

IF 

HKS 

T.rg.4) 

.AND. ( J.GT. 

1 1 3 )G0T0342 

' 4 ? 

IF 

( f.K’:. 

T.fiO.S) 

.AND. ( J.FQ. 

113 )G0T02 

f {') 

IF 

1 ’ 

I { K S 

T ■■ j 1 

r,F>.i.5 3 

7 

» A ’^1 0 • C J ® (ji T ♦ 

11) )GPTn34? 

-1 '' s f 

r;! 

:• CA 

1 ! j t K, 

1, L T I 

CK(M,MA 

f F , F 1 3 


i"'} 

F-O 

T-f.cn 

T=i ,N6 



(ii '2 

3 • ') A 

FP TF 

T}= -F 

1{NT3+S{MT3 


e> *'"/ 

C.A 

!..L S 

UR{.AFP,! 

MA,H.ZA3 


t b 

nil 

3'}lh 

Tsl ,NA-. 




A F p ( M ) = F 1 ( NT ) « U .-F 1 ( NT ) ) 
CALLSU13{AFPrNA,HtZ0) 
r_Ai,isua{s,NAfH,zp) 

ZH = ZP-'.7.S. 

sa»=l ./;{ 

ZH= (-25.*F| ( 1)+48.’»F1 t2>-36 
ZR* ( -23 {, 1 ) +48.,’S*S 12 3-36 


,NAfH,ZP) 

FI ( I ) +48.’»t-l (2 > -36.^F 1 ( 3 ) + 16 .^F 1 ( 4 )-3,=«'F 1 rS ) J'S'HH 
S {1)+48,’S*S 123-36. S'S {33 + 16. *S (4)-3.4S { 5 3 3NeHH 


n 

3 02 

, If 
If 

{ FvST.{4t .,4,3 ZQ; 
tF'5T.'^NE,4 3ZM:' 

;r', • 

■ ni/i 

3 33 I M 0=1:., 5 

"''0( 


iMTg,p,f>f>,:sw 


■R^i 

!0T4,Zn ■ 

j 

PR 

INTT.ZR 

3] ! ■ . 

Rh; 

I NT 6, Z A 

^31-? ^ 


ior/,zo 

013 ' V, 


! NT .'J » / P 


G004 


I Sf^i 


S '0 U ,R C t: z T T if- ^ J 

3 1 4 


If- (0 ST, !,'-”„/( ) pnifvIJRq 

1 1 7 


n»r.Tio 



EFA = 0,. 

? 1 


'D(i304N=l,!\j/\,To 

3 27 


*"Hrrn,i,FTA,nN),Fi: 

,, 

3 0 

A f TA=rrA+.5oocoo 

*'5 7 3 

3 3 

COOT I our \ 

3 27 


GOTOiCT. ' 

3 3/3* 

1 0 5 

GO TO 3 31 ■ ■ 

"5 3 I 

33 1 

cr:-MTiNur ' ... 

:i 

1 /, 7 

IF ( ^' ST. F0.7 ) GCT0321 

'•> 3 c.^ 


I i- { K S T . R ) G0T03A1 

3 ] 


If ('“'ST. ::t,..A)GGTC;3 51 

4 /.. 


r r ( ,7 '• T . ri)^ q j r, 070301 

■'. 4 



0 3 


7 

.'■ '"5 !, 


I *■■' 7 = 3 

' ^ 4'* 


:-'i:rri/,7 

vr3 

? ] 


i > 'Z^'* 


7 ::i T = 3 

5 b 5 


!’ = y . 

", S ( • 


r K 

■> *■> 


roiTriA'T 

iFTJ 

,j M 

!<■ T = ■ 

■1 6 1 


F. >/ = . 1 

BO? 


P 1.', = . ■? 

•i (-1 ■’■ 


10^=3 

■i <') A 


nr. ■^01 *'7' 

■: h B 

'3 h I 

'< 1. T = A 

f} 


!.B = . 1 

,' 6 Y 


■■»B = ,U3 ■ 

■r/i; 


I K = '3 

■1 n 


GC^OIA? 

.7? 

34 1 , 

STi/P 

> 7 ? 


r,F-’0 


FORTRAN SOURCE LIST 


MLC.nO'' ' ■ 

: NO FOR above' ASSFI^BLY 


IBMAP ASSEf-^BLY 


// 


FQPTRAM SnURCt LIST 


STt'.TF^''^^T 


'•EGaO^ 


S i I (' ;■ n ! : T p.! 2 F L I C K { i-! » M A » AF , D 1 F 1 0 2 F ) 

AF(21D) , C 1 F ( 2 10 ) t G?F ( 2 10 ) 

CuLLFLiF'-MBtO) 

A 4 . ]_ 

Hl. = l . /( I2.«H) 

Flil = l ./( 12-*{H**2) ) 

H.j 1 = I . / ( 1 2 . ( ( ? . ) ) 

HK = ( U/R. 


rn 4^ L-l , j. ftc f I 4.7 \ +T * A,F ( L+R I~3.^=AF ( L +4 ) ) *HF 

1 4 f ' 1 « ■=^" A f * ( L 1 ^ ® ^ ^ ^ ^ ^ ^ A / % 

Y ? = ! i . n T ( 4 5 . ='^ A F { L ) - 1 5 4 . A F { 1. + 2 ) + ? 1 4 . « ^ M L + 

1 -iF6.*AF{L + 6)+6i.4AF(L + 3)-l-*''-A. IL . ■<) ) 

'i::>2F (i )=(HK*Y?-Yi )/HKl 

y|A!/o;«p’F a-?)-8.=!‘AF (L-l I+a.^t'AFIL + l )-’• .*AFa+2) ) 

DiF d )=(HK--!'X2-Xl)/hKl + ) 

( [ ) = (HK*Y2-Y1 5 /HKl 

DO'«:l.>‘ir<,'<A ..C11-1H.3A »Ai:lL-?|-16.*AF(L-3H'3.«AFa-4) I 

*F.:HT*r;‘0«HU)-;e;«r(L;2)*36;«Fa-AI-l6.«F<L-6,*3.*AF(L-») 

l;i:i^Hrir«!^rFa!-13i.«FU-l).21A.*AFlL-2.-l«.^AFa-31 

_?j+? 14 .’!fAF(L- 4 )- 156 .*AF(L- 6 ) 

Y/»lt j 1*^' t45»*AF ( L ) -l'>4.*Ai ( L 2 ) - ■ 

1 4-6l.’*AF{L-8I-10.=5tAF (L-iO) ) 

. : 'U2F LI- ^ ■ 

'-'rETURL :4 ^ : , //'V'v ' 


I8MAP ASSEMBLY FLICK 


MO M!'SSAi.?t:"3 


F fi R A F 0 V F ASSEMBLY 






73 


3T;\Tt-iVTr!T 


FORTRAN! SnijRCE LIST 


0 

$ IB FT 

C ;< c T T 

1 


F!J^''RruTr-;r KUTTA (NA, P'', 

2 


lC5IC:iL 1. •NG 

3 


■■ 1 >5::NS! CCAF ( 2 IC ) , S( 210 



li.l S (2 10 ) , U?S ( 210) ,Av'( 3 



■- » D.K(6),CL! 

3 

/„ p 

FljR>''A-(4315.5) 

CALL FLlJlOOCn) 

(■. 


CAU.FLiCK (H,N1 A,AFtDIF,I 

1 


C.ALLf iC!<(H»l:A,S,0iS) 

IT 


no: 1.1 L=i,NjA 

11 


n3F{L}=-(i,F{L)*D2F{L)+[ 

I 2 

1 1 , 

02'^ (L )=-PR*{AF{L)«DlS(!. 

1^) 


ni' ( u = 0. 



■' </{ 1)=1 . 

1 6 


AVVt 1 )=1. 

IT 


Ay!2)=i.2 

il, '■ ■ 


A V V (?) = !,. 

21 


AV( 3)=1. 

? 2 


AVV(3 )=1. 3 

2 3 


AV(4 )=1.,2 

2 A 


AVV{4} = l..;i 

25 

2 6 


■H' { 5 ) = 1 . 6 

Ayv(5)=1.2 

2 7 


Ay (6 ) = !.,« 

3 0 


ftVV{6)=1.2 

3 1 


.AV(?)*1^6 

32 


AVy{7) = l.5 

3 3 


AV{3!=2. 

34 


.Ayv{2)=?.i 

35 


NS = 3 

36 


HLlUsO. 

3? 


FiF ( j. ) ■=(.! ft 

AC 


CPHDsQ. 

41 


COM 1 )=0. 

42 ^ 


.11^=1 , ■ 

;,4,3 


MK-l ■ ■ . , , 

,44/''' 



45 


K*,i;.O::'\ 0'' / ■ ■■ ■ 

46 



47 

’V- 

■i,= l 

50 



51 


T (11 = G„ 

52 


R=lv»-r 

5 3, 

XT / 


p = 0. 

*t 

3T 

^ 'fy 

i.., L,! 4 4 4 I 2 f 5 

56 


CALL FLUn(lOQ) 

57 


A'=0. 

6 -r.' 

6 


IF n .G7.2)A=1./2. 
l!'-( I .CT.4)A=1. 

L 6 


A1=!!11MLH-A=5^D2F(L1*H 

/■7 


A,'-r'2F (1, )-t-A*C3FtL)*H 


> » I i f I 


NH:r:,l'r=(L ) A * C 2 S ( L 1 « H 

AA-. “I- (L 5 + rilf" (!,) *H«A 



FORTRi^f'' SfiURCC LIST KUTTA 


j >r, STATS MSfJT 


7 ? s ,; = S(1„ )+,■)! S{L)*H* A 

7 3 ^'=<>{1 ) + !:'K { 

7S ¥ = •> + ■ LI 1-1 )*A 

73 7 

7 6 i; = T f L l + DM ( r -U*A 

7 1 V = ;-M-rr ( i-i )*A 

ICi,: rA(nxH«'-Y 

l,Cl iJL(n=i!*Z 

ino C-'( t 5 =-7^:;- ( X7/A7+fiA=5--7 + C« (U-?.^A1*Y ) + ( A5-A1**? ) i 

!03 

5 0 A - n ( 1) =-H* p R * ( A4 * V A 3=^^ X ) 

ir,'i AAA CrMIiSU;- 

’ ST ■■ ( L + 1 )=C(L) + l./6.*(i}K(2)+?.*DK(3)+2.=^'-;jK!45+DK(5) ) 

! 1C P = P-i-i./6.X={DL!?)+?.*DL{3)+2.*i.L(4)+DL(S) ) 

1 I 1 C = . + 1 . /6 . ’M C^' ( 2 ) +2 . *D*' ( 3 ) +2 . *0!^' ( 4 ) +D?' ( 5 ) ) 

I 1,7 T { L>n = T ( L) + l ./6.* ( D^i{2)+?.*^r:( 3)+2 . YV:;rv(4 >+DN( 5 ) ) 

r 1 3 = ii 4- 1 . / 6 . ’M DU { 2 ) + 2 . *00 { 3 ) + 2 . *DO ( 4 ) +DD ( 5 ) ) 

114 Tl = T{L+n 


i, IS 
5 14 

! I Y 31 
I 20 

L21 32 

12 ? 

I 2 3 
1. 26 
12 / 

I 3i'l 

i31 36 

1.32 
13 3 
1 34 

1,36 37 

1 4 0 , 

141 

142 ■ , . 

1, 4 3„ ■ 
■144c 

' :x ■■■Y; : 

^"■"c 

'■1,47' 

130 
1 61 
1 52 

I C 3 
1 

} 5 5 3 5 

1 56 
7)4 


u * j - '-.i ( L + 1 ) 

IF{L-UA )31,32,3? 

L = L + 1 
GCTU33 
f T T ( '3 A ) 

LAG = AH5 { P ) . LT. .00001,. AKD.TT .LT . .00001 

IF(tAG)G0Tn35 

AT(!I)=TT 

Ap{ 1 ! j-p : 

IF ( n-NS) 36 t3T»3 7 
I I = l + I I 
W = A.V(!I) 

WW*.AVy(in; 

■G[VTn3C ' 

IF {KK .tti.2)30T035 

r, 1 a ! AT ( 2) -AT ( 1) ) / ( AV{2)-AV(i)) 

C2=fAP(2)-AP( 1))/ lAV{2)-AVtl) ) 

CC^ = ( AT( 3) - ATI 1 n / ( AVV (3 )-AVV I 1 ) ) 

C4 = ( A P { 31 - A P 11 ) ) / t A V V I 3 ) - A Y 1/ { 1 n 
0FP-=- A,P '(1)c ,' ' ' ■ 

(1 )C. '■ ■ ■ 

^XC=,Ci4C4-C2*C3. 

0 A V « I C4 *5* OF T-DFP*C3)/XC 
OAVV=IC14riFP-C?*DFT)/XC 


WsAVlD+DAV 


V,.w = AVV( 11+DAVV 


KK=YK + 1 


GOT 03 3 


return ' 



IBMAP ASSFY^SLY KUTTA 


0 FfCSAGl 


roi 


Al'.nVF ASSEPBLY 



[ 


‘•tG 


* 'll ■''' 


p 9 n 


cno 4 
T SM 




0 f . n-tFTCGlin 

1 

7 


3 

4 
3 
6 

7 

10 

12 

1 

I 5 
1.6 
] 7 


,i;V SUO{Y,M,H,AREA) 

fJiMFNSlOM Y(210) 

S 0 ;-'~V = 0 . . 

C i f/ ■ r .-v 


K = ^ I -1 : 

004 l = 2,.M,''2' 

n , ifEV = Gti : Rrv -- s - v .; r ) 

nC'jI = 3 , n ,»2 

5U‘0:l;=SUYG 0+Y( I ) 
^■'=!i/3.*(Y( I )+4=>-.Sl;- Gv + 2 


0 " Tij ;7 

F ! : ?") 


J \ 


FORTRAN SOURCe LIST 


. " SUOuil+Y { f ' j ) 

I SWAP AS SEP FLY UD. 


R'r; Rf^SSAGC^ FOR above assfnbly 

OA ' 

IS BEING r^TERFD INTO STORAGE. 


IBLflR — JOB OCOOOO 



7 ^ 




iiPPENDIX 5 



77 


E6004 


SCURI 


SlATEi^FNT 


FORTRAN SOURCE 


01 O!- NS I ON ST{ 15) ,FN{ 10) , JCI (25 ) , JCF (25) » JQ( 25 ) ,C( 2 j ) , 

I AF ( 100,25 ), 0( 160,25 ) f AETA { 160) ,JEQ{ 25 ) , 

'c- NNV( 10) ,SST (15,5 ) , FFN ( 1 0 , 6 ) , J JC I (25,6) , JJCF{25,6) , 

3JJh.-H25,6) , LCH(6) ,NIB{6 ) ,MFo{6) 

LOGICAL LG 

FfjRNAT ( 5 X, 30HS0LUT ION OP-TA I NFi) 0 ITH PARA . D I FF . OP THOO ) 
F0RMAT(2X,6HAF{3)=,E13.6,6HAF(6)=,E14.4) 

FORMAT(5X,45HSOLUTION OBTAINED DIRECTLY FRON ORIGINAL -ON.) 
FnRYAT(5X, 10HPARAMET>R= ,Fli.4) 


7 

100 

Ft,JRMAT(F4.2 ,F4.2t2I 3) 

10 

2-) 9 

FOSMATdl) 

11 

9U4 

F0RFAT(fc!Ell,4) 

12 

2 3 

rORMAT(4HKHP=, 13) 

13 

14 

FORMAT {2F 5, 3) 

14 

15 

FORMAT! 711) 

15 

l6 

FORMAT ( 31 1 ) 

16 

13 

F0RMAT(3F5.3) 

17 

11 

FORMAT (3 12) 

20 

10 

FORM AT (2 I 3, F5. 3 ) 

21 

117 

Fl}RHAr{F5.3,13) 

22 


RGADl 17,DB, JA 

24 


RGAD1C:,NKHP,NR,H 

27 


REAOll, (Nrw( I ) ,I = 1,NKMP) 

34 


R6A01 1, {NIB( I ) , I=1,NKMP) 

41 


RE ADI 1 , ( MFB ( I ) , I = 1 ,NKMP ) 

45 


nU121=L,NKMP 

47 


M^NIPII) 

50 


MN = MFiMn 

51 


N:=NNV(I) 

52 


R6AD13, (SST(K,I ),K=1,M) 

57 


REA014, (FFMIK, I ) ,K=1,MN) 

64 


READIS, ( JJCI CK, I ) ,K=1 ,N) 

n 


READ15, ( JJCf (Kfl ) ,K=1 ,N) 

75 


READ! 5 , ( J JEQ ( K , 1 ) , K = 1 , N ) 

103 

'/12V.: 

. CONTINUE- : 

105 


READ 1 6 , { i CH ( K ) fK~ It mm ) 

112 


mmzo9,m ■ 

1 14 


AETA{1)=0. 

115 


NAA=NR-1 

116 


D04A4k=1,NAA 

1 17 

/.| 4 4 

AETA{K+1)=A£TA{K)+H 

121 


nr;!i7KMP=i ,nkmp 

122 


LG=( (LCH(KMP) .Eq.2) .OR. (L 

123 


NV=NMV( KMP) 

124 


IOC=NiD(KMP) 

125 


IFE=NFB(KHP) 

126 


D013J=1,HV 

127 


JCI ( J ) = JJCI ( J,KHP) 

1 30 


j[:u( j)=jjr-o( j,kmp) 

131 

IF 

.JCF( J) = JJCF{ JtKMP) 

1 3 3 


0019.1=1 ,11'!. 

1,34 

19 

ST ( J)=SST( J,KMP) 

135 


nO?uOJ=l, IFB 

1.37 

200 

F;;( J)=FFM( J,KMP) 


7 ^ 


GO 04 
ISM 

141 

142 

143 
146 

151 

154 

157 

162 

165 

170 

171 

172 
1 73 
I 74 
175 


201 
2 02 
2C3 
2 04 
205 
214 
22 3 

2 2*5 
2 26 
2 30 
2 31 
004 


SOURCE STA^EMFMT 


FORTRAN SOURCE LIST 


CALL ANT( NV» JEQ» JO , IC ,MI , JC , JCF ) 

LLM=Lr,H{KHP ) 

IF{LCH{KHP).E0.0)CALLBST(AF,NR,H) 

2K|if>tM'Ic!'c!jc?LG!LLM^ ’ f*’ ' NR. S T, FN, JCF , JEO , JC I , JQ, i c , M I , 

I f- 1 LCH{ KMF ) . E0,3 ) GOT020 
I F ( LCH{ KMP ) .EQ.O ) G0T0i7 
IFILCHIKMP) .EQ.UG0T021 

J LCH(KMP) .EQ.2 3G0T022 
22 PRINT302 

PR IMT23,KHP ' 

GnTrj304 4' 

2 0 0iJ24MS=l ,Np, 

1. tj 4M J — 1 y \J 
0 ( MS tNJ ) =AF { N'S, MJ ? 

1 » 1 G^L LM J ^ ^ ^ ‘ ^ ^ » I I » R MP, AF , Q, C , JC 

GO TO 3 04 • 

PRINT 301 

Pk I f'-’T 3 JjC f B 

00805 N=1tNV 

rFIKpn'l^rr. o. f APfl ,N) ,AFTAC I),1 = 1,NR) - 

COMTIMUE PRIMTg04, (AF(I,M) tAETA ( 1) , I = i ,nr ) 

GO TO 17 
CONTINUE 
STOP 
ENf) 

IBMAP ASSEMBLY 


24 

21 

304 

B 0 5 
I 7 


NO YF:-»SAGtS FUR ABOVE ASSEMBLY 


7^ 


004 

ISM 


FORTRAN SOURCE LIST 


SOURCE STATEUfi'iT 


0 stbftc mist , 

1 SU"^ri:T!NeMIST ( JA,DB,H,B,NV,NR,ST,FN,JCF, JFW, JCI , IC,MI ,KMP,AF 

1 ,Q,C, JC,LG) 

2 LOGICAL LG 

3 or Cr, SIGN JCF ( 25 ) , JEOt 25 ) , JC I ( 2 5 ) , JQ ( 25 ) , J C C 6 ) , AF ( 160 , 2 5 ) , 0 { 1 

260, 2 5), 01(160,25) , AFF ( 160 ) , OIF ( 160 ) ,C(20) , AF 1 { 160, 25 ) , ST ( 2 5 ) , FN ( 
310) 

4 rn9r;9i = i,ic 

5 909 C{ I ) = 1.CT<0D 

7 [•:VV = 'V-1 


10 


rc^ 53 2 1 = 1, MW 

11 


TF( J&{ I ).E&.l)GOTG512 

14 


D05r;;N=i,r4R 

1. b 

513 

AFF (N ) = AF (N, I ) 

17 


CALL TICK (H, MR, AFF, OIF) 

20 


D0514N=1,MR 

21 

514 

AF (M, I+i )=D1F{N) ' 

2 3 

512 

CONTINUE 

2 5 


XYZ=DB/6. 

26 


D G 8 = (i B / 2 • 

27 


b=0. 

30 


JJ = 1 

31 

H07 

DO 80 1 1 = 1, NR 

32 


00801 N=1 ,MV 

3 3 

801 

AFl ( I ,N)=AF{ I ,N) 

36 


D=B 

37 


AA = . 1 

40 


CALL RKM ( H , 0 , MV , NR , ST , FM , JCF , J FQ , JC 1 , JO , I C 
1,LG) 

41 


D0802H=1,NV 

42 


D0802I=1,NR 

43 


CALL FLUNOOOO) 

4 4 


AFl ( 1 ,M) = AF ( I ,N)-t-Q( 1 ,N)*UBB 

45 

B0,2 

CONTIMUE 

•50 


D=B+D8B 

51 


CALL RKMIH, n,NV,NR, ST , FM , JCF , JEQ , JCI , JO, IC 
1 , LG ) 

'■5t 


D0803N=1,MV, ■ 

^^53 


008031=1, NR 

■: 54 :. 


CALL FLUN(30Q0) 

;55-'' 


AFl ( I ,N)=AF1( I ,N)+Q1( I ,N)*OBB 

56 


0(1 ,M)=0{I ,N)-»-4.*Ql( I,N) 

57 

803 

CONTINUE ' 

62 


D=B+DO 

63 


C A 1. L RKM { H , 0 , N V , MR , ST , FN , JC F , J EQ , JC I T JQ , I C 
ItLG) 

64 


D0B04 N=1,NV 

6 5 . 


00804 1=1 ,MR 

66 


Q(I,N)=0{ I,N)+Q1(I,N) 

67 


AF( I ,N)=AFt I ,N)+XYZ=4'QC i,N) 

70 

804 

CONTIMUE , ' '■ ■ , 

73 

904 

F0R«AT{6E15.5) . . . ’ 

74 


B=e+n8 ' , 

75 


JJ = JJ-H 

76 


IF { JJ.EC. JA )G0T0808’ 


FORTRAN SOURCE LIST MIST 


004 

ISN 

101 

1D2 808 

1C3 
004 


NO MESSAGES 


SnURCF STATE'-’FNT 

GCTOaO? 

RETURN 
END . 

I8MAP ASSEMBLY HIST 


FOR ABOVE ASSEMBLY. 





;oo4 


FORTRAN SOURCE LIST 

ISN 


SOURCE STME;-’F?;T 

0 

$IRFTC HIX 

1 


SUPRPUTIME HIXC H ,0, NV,NR, ST,FN, vJCF , J FQ, JC I , JQ , I C , N I , KMP , AF, 0 , C , J 



1C,1G,LM) 

2 


LOGICAL LG 

3 


r: ! :- F N S I ON S T { 2 6 ) , JCF ( 2 5 ) , JEQ I 25 ) , JC H 2 5 ) , JO I 25 ) , JC { 6 ) , AF { 1 60 1 2 5 ) , 



10( 160,25) ,FN( 10) ,C(20 ) 

4 


IF{LH.e0.2)G0T03 

7 


G0T04 

10 

3 

NN=1. : ' 

11 


D01I=1,NV 

12 


IF{JC1{I),FC'.1)G0T02 

15 


GQTOl 

16 

2 

C(NKt)=AF{l, I ) 

17 


NN=NN+1 

20 

1 

CONTINUE 

22 

4 

AA = .Cfl 

23 


CALL RKNt H,6,KfV,NR,ST,FN, JCF, JE9, JCI ,Ja, IC,NI , K M P , 0 , AF , C , J C , AA 


' 

1,LG) 

24 


RETURN 

25 


END 

004 


^ IBMAP ASSEMBLY HIX 

NO NES 

.sage; 

> FDR ABOVE ASSEMBLY 


FORTRAN SOURCE LIST 


‘^■ 2 - 


GOOA 
. I SO 


SOURCE S' 


AT:P';-;r';T 


0 SIBFTC ANT 


1 


SUBROUTINE ANT ( NV, JEQ » JQ» IC I , JC , JCF > 

2 


DI‘''HNS!L!JJE0(25) , JQ ( 25 5 , JC{ 6) , JCF { 25 ) 

3 


o 

m 

il 

4 

8 

jv. (n=JFv.un 

6 



7 


PGY I=lr’-'S 

10 


N=I + T 

1 1 


IF{ JEO( I) .EC.l) JEQCN)=2 

lA 

7 

CONTINUE 

i h 


00117N=1,5 

17 

11.7 

JC(N)=0 

?1 


IF-(MI .EU.0JG0Tnil6 

2A 


D01IN=1,HI 


11 

j(:{N) = i 

?7 

1 1 6 

I A=n 

30 


DD9I =1,NV 

31 


IF(JCF(I).E0.0)C0TG9 

3A 


IA=I A+i 

35 

9 

CONTINUE 

37 


IC = IA 

AO 


RETURN 

A1 


hnd 

sonoA 


IBMAP ASSEMBLY ANT 


MO HESSAGES FOR AfiOVE ASSEMBLY 


t 


FORTRAN SOURCE LIST 




G004 

ISN 


SOURCE STATFMEf-.T 


0 1. IB FTC RKM 


1 


SUBROUTI ■ :FRKi';(H,B,NV,NR»StFN, JCF, JEQ, JC I t JQ » I C , M I , KHP , Or F , C » JC 


1 

,L0) 

2 


LilGir.lL L3 

3 ^ 


LuGiC;.-.L LAG,LAn(6) , , , » 

4 


l)^-H-?i:linri5(25) fFNllO) , JCF( 25 ) » JEQ t 25 ) t JCI ( 25 ) , JQ I 25 ) t F t 160.251 


T 

x 

0{ 160,25) .;iQ(5) ,0t Ilf 11) fC{20) ,G{ 11 5 ,DFI5,25) ,FFt 160) ,DFF ( 160) 


i... 

YI Il.lli fldl) ,LA(ll»ll)fCYIll»ll)fX( ll)»JC(6) fLSUMCll) 

5 

200 

F0RYAT{E14.4) 

6 

212 

FnRMAT(2Xf 4HG( I ) ,E14.4) 

7 

?ll 

FnRlCU(4r3) 

1 0 

213 

f URH AT ( 2X , I 3 f 7HF ( 1 , I Z ) , 2 1 3 ) 

11 

214 

FO RM A T ( 2X , 2 HD F , E 1 4 . 4 , 2 I 3 ) 

12 

21 5 

FORM A T { 2X , 8H Y { J , NN ) = f E 1 4 . 4 , 4 1 3 ) 

13 

216 

F0RYAT(2 X,9HLA{ I ? , IN ) . E 14. 4 , 2 I 3 ) 

1> 

217 

FnRMAT(2X ,8HLSUN( IP) t2I3) 

15 

21 S 

FnRYAT(2X. 8HCY(NtLL) ,E14.4,2I3) 

16 

219 

FORMA T(2X,8HCY( I ♦ JK ) , E14.4, 213 ) 

17 

220 

FORMAT (2X,4HXl I ) ,2E14.4) 

20 

221 

F-0RMAT{2Xf4HC( I ) .2E14.4 ) 

21 

222 

FGRMAT(2X,2HI=, 13) 

22 

400 

FORMAT! 1X,2HDF,7E14.4) 

23 

33 

FORMAT (2F5. 2) 

24 

32 

FORMAT! 3F5. 2) 

25 

31 

F0RMAT(10X,7I2) 

2 b 


CALLFLUN! 3000) 

27 


PRtNT31 T ! JEU! I ) * i=l»MV) 

34 ' 


PRINT31f {JCF{!),I=1.NV) 

41 


PHINT31 , ( JCI ( I ) » 1=1 »NV) 

46 


PR1NT31,! JQ!n»I = l.NV) 

53 


PRINT32,S!1) ,S!2) fS(3) 

54 


PRIHT33.FN! 1 ) fFN(2) 

55 


PRINT200,A . 

56 


PRINT211»KMPtNVf ICf HI 

57 


B0.!l)=-1 . 

60' 


80{.2I*.0'01 

6l 


■8Q{3) = .*0001 



B'0(4l».Q0001 , 



8Q{ 5)=.OOC001 

64 


PRINT221*C! 1 ) tC(2) 

65 


ISUH=0 

66 

461 

Z=A*2. 

67 


KK = 1 

TO 


KKK = 1 

71 


IF ( IC.EG.0)G0TC4ia 

74 


0010 I=lf ic 

75 


DOlO J=l. IC 

76 


IF! I .E0.J)DII.J)=Z 

101 


IF! I.NE.J)0! IiJ)=0. 

104 

10 

^ CONTINUE : 

10 7 

418 

JK=1C+1 

110 


KK = l 

1 u 


KKK^l 

1 12 


ZU=1 .COCOO/6. 00000 

113 

100 

00320 J=l. JK 


FORTRAN SOURCE LIST RKM 




G004 


ISM 


114 


115 


120 


121 


124 


127 

11 

131 


132 

4 0 

133 


134 


135 


140 


143 


14A 


147 


150 

13 

151 


15? 


1 5 3 

15 

1 54 

12 

156 


157 


160 

19 

162 


16 3 

547 

165 

21 

166 


167 


170 


171 


174 


177 


202 


205 


210 


213 


2lfe: 


21-7^^'. 

501 

2 20' ^ 


221- 


222 

20 


2?5 

?2h 

221 


2.32 

235 739 

237 

240 

243 

244 300 
2 4 

250 

251 
2 54 


sourcf: statement 


C MLFl.UfH 3Cr;0) 
ir( IC.EC.COGCTGAC 
DOll 1 = 1, IC 
IF{J*E0.1)G(I)=C1I} 

IP{ J.GT,1)G( n=C{ n+D( J-ltl) 
CONTI '‘IJF 

PSINT221 ,GI 1 ) ,GI2) 


nci2 IZ=1,NV 

IF( JO(IZ) .E0.1)GGTC15 

IF(JCI(IZ).EQ.O)GOTC13 

IF{ JCH IZ ).EG*1)F( It IZ)=G(N8) 

NB=NR +1 

007012 

F{ 1 , IZ)=S{MS) 

NS=NS+1 

G0T012 

F ( 1 , 1 Z ) =0 . 

CONTINUE 


K =1 

0019 M=ltNV 
DF( 1,M>=0 
DO 54 7 M= 1,5 
0F{M,1)=D. 

D020 M=2,5 
N 0 = 0 

0020 MA=2,NV 
CALLFLUN{50D) 
IF(N.EQ.2)A=0. 


IF( M.FQ. 3 )A= 0.50000 
IF(H.EQ. 4 )A =. 50000 

IF (M.E 0 . 5 )A=i . 00000 ^ 

IF { JFOI ha ) . EQ.O) DF ( H, MA) = H* (DF { M -1 ,HA +1 KA+F I K , MA) ) 

I F ( JFQ( HA ) • EQ« 1 ) GOT050I 

IF( JFaiMAI.EO.Z) DF(H,MA)=DF( 1 ,MA) 


,G 0 TG 20 ■ . , 

'NO=N O + 1 ' 

CALL FUNSON(H,Bf A,KtNV,JOtH,MAtQ,FtDF,DFPtKMP,NO) 
DF{H,HA)=DFP 
CONTINUE 
.D 0739 H= 1 ,MV 

If\‘ JO wlE0^0)F(K+l,M)=F(KtM) + ZU*(DF{2,M + i)+2.«DF{3tM + l) 
1 + 2 .*DF( 4 , H+ 1 ) +DF ( 5 t M+ 1 ) ) 

If ( JC(M).e 0 . 1 )F(K-fltM)= 0 . 


CONTINUE 

-f* 2 . 

IF{K!eQ.NR)GOT0300 


G 0 T 021 
NN = 1 

IF(IC.EO.O)GOT 0410 


0023 H= 1 ,NV 

IF ( JCF(M) .E 0 - 0 )G 0 TQ 23 

IF ( JCF( H) .EO. 1 )Y( JtNN)=F{NR,H) 



{ iEGn04 


I SN 


SOURCE 


tr-^ENT 


SOURCE LIST RKM 


257 
2 60 
• 261 

263 

264 

2 65 
2 66 
26 7 
272 
275 

300 

301 
.3 02 

3 03 
3 04 
305 

310 

311 
313 
316 
320 
3? I 
324 
327 
332 
33,5 
340 
34 3 
346 

351 

352 

3 53 

354 

355 

356 

36 1 

362 
3 ^ 3 :, 
364 

'v'.3,6'6 ^ 
: "367 
374 

401 

402 

4 0 3 
405 

412 

413 

414 

415 
' 416 

421 
4 22 
42 3 
425 


PC I NT 21 5 5 Y ( J , , j js, 

0 i',' + 1 

2 3 CniiTINUt: 

'.■:4:'!l IP=1,5 

3'.7(D1 IN=i,IC ' i 

T ( I'i)=ABS (¥( J,IN)-FNI IN) ) 

LaG = T (IK!) ,Lt,8Q( IP) 

IFILAG) LA(IP,IN)=1 
IF(.N0T.LAG)LA(IP,IN)=0 
4ul CONTINt..iP 

1)0402 IP=1,5 
ISUM =0 

00402 I!j=I,IC 
ISUiM = LA{ IP, IM) + ISUM 
LSUM{ IP) = ISUf*' 

402 CONTINUE 
D0403 IP=1,5 

403 LA!){ I P)=LSUM( IP) .EQ. IC 
IF-(LAl}( 5) >0010410 

320 CONTINUE 

LAG=KK.E 0.2 

IF ( (.NOT. LAD(1 ) ) .AMD. (LAG). AND. I JC( 1 ) .EO. 1 ) )G0T0460 
IF C ( .N0T.LAD(2) ) .AND. (LAG). AND. { JC { 2 ) .EO. 1 > )G0Tn460 
IF C {.NOT.LADO) ) . AND. (LAG). AMO. ( JC ( 3 ) .EQ. I ) ) GDT0462 
IF{ (.N0T.LAD(4) ).AND. ( LAG ) . AND. ( JC ( 4 ) .EQ.l) )GnTn462 
IF (LA0(5) .AND. ( JC(5 ) .EO.l) )GnT0410 

IFC {.N0T.LAD{5) ).ANO. ( LAG ) . AMD. ( JC ( 5 ) .FQ.I) )GDT0463 
IF(LA0(4) .AMD. ( JC(4) .EQ.l) )GOT0410 
IF {KK.6Q.2)GOT0410 

0025 L= 2 ,JK 
D 025 M= 1 ,IC 
IL=L -1 

CY(M,LL)=(Y{L,M)-Y{1,M) )/Z 
PR1NT218,CY(M,LL) ,M,LL 

25 CONTINUE 

0026 1 = 1 , IC 

CYdf JK) = FN(I>-Y(l,r) 

PRINT219,CY( I , JK) , I ,JK 

26 CONTINUE 

CAL'L AiEO( JK,CY,X) 

PRINT 220 , (X( n ,1 = 1 , IC ) 

PRINT 221 , (C(n,I = l,IC) 

D 070 I= 1 ,IC 

c( i )=c( n+x{ I) 

70 CONTINUE 

PRIMT 221 , (C{n,I = l,IC) 

'kk=kk+i , . 

JK = 1 

GOTOICO ■ 

410 0041=2, NV 

IFIJQID.EO.I) GDT05 
G 0 T 04 

5 D 06 N= 1 ,NR , 

6 FF {N)=F(N, I-l ) 

CALL TICK{H,NR,FF,DFF) 



; ' ' ' ' ' ' 

£G004 , ^ , 

1 ISN 


SOURCE STATEMENT 

1 


D07N=1,NR 

427 


F(N, I )=DFF(N) 

4 20 

7 

CfH'iTIKUE 

432 

4 

COMTIAiUE 

434 


GOT 04 79 

435 

460 

IF(LG)A=.Q05 

■;| 440 


A=0.025 

!| 441 


0010461 

' 442 

4 62 

ir(LG)A=.002 

• ‘ 4 4 5 


A=.005 

- 'f 4 6 


GCjT0461 

447 

46 3 

IF (LG)A=.001 

; 4 5? 


A=.00? 

45 3 


GQT04fcl 

■ ^54 

479 

RETURN 

: 45 5 


END 

tG004 

NO MES 

SAGES 

FOR ABOVE assembly 


FORTRAN SOURCE LIST RK.f^ 


IBf-lAP ASSc^'BLY RKM 







FORTRAN SOURCE LIST 

I I SN SQURCF statement- 

I 0 SIOFTC TICK 

i 1 SUBROUTINE T I CK ( H , NA , AF , D1 F > 

? OIMFMSION AF(?20) ,D1F(210) 

3' CALLFLUNI 3000} 

A NHA=NA-A 

5 NB=NNA+1 

6 HH=1 ,/( 12,«H} . 

7 HJ = l./{ 

10 HK=n../2. }**4 

11 HKl=HK-l. 

12 D046L=1,4 

13 Xl = (-25.=<'AF{L)+48,=f'AF(L+l }-36. 4AF { L+? ) +16. «AF { L+3 }-3. *AF (H-A) )*HH 

14 X2 = HJ*{-2 5.*AF(L)+48.*AF(L'+2)-36.*AF( L+4)+lA.*AF{ L+6)-3.*AF(L+C) } 

15 46 OlFIl ) = «HK=f^X2-Xi)/HKl 

IT DD47L=5,MMA 

20 Xl = ilF4(AF(L-2}-8.*AF(L-l)+6.’?'AF(L + l )-l.^AF{L+2) } 

. 21 X2 = HJ4-{AF(L-4)-8.=«'AF{L-2)+8.^AF(L+2 )-i.=«'AF(L+A-) ) 

22 47 D1F{L)={HK*X2-X1)/HK1 

24 DD4P.L=NP,NA 

25 X1=HH«{25.*AF{ L)-48.*AF(L-1 )+36,*AF(L-2) -16.*AF(L-3)+3.«AF(L-4 ) ) 

.i 26 X?=HJ«(25.=«=AF(L)-48.=S‘AF(L-2)+36.=!'AF (L-4)-16.=?'AF(L-6)+3.*AF(L-F) ) 

I 27 48 niF(Ly=(HK*X2-Xl)/HKl 

i 31 RETURN 

■i? 32 END 

;|EGn04 IBMAP assembly TICK 


I NO MESSAGES FOR ABOVF ASSEMBLY 






1 MO Mt-:SSAGeS FOR AOOV6 ASSEMBLY 

I 


1=0004 
. ISN 


SOURCE STATEMENT 


FORTRAN S3URC' LI 


N'EGOf 


$IBFTC BST 

SUBROUTINE BST(AFtNR,H) 

DIMENSICN AF{ 160,25) 

ETA=0. 

002211=1, NR 

AF( I ,1)=ETA + EXP(-ETA)-1. 

AFU ,5)=1. + {1.5-1. )*EXP(-ETA) 

221 ETA=nTA+H 
RETURN 
END 

IBMAP ASSE'^RLY 


NO messages for above assembly 

MFC no 4 


IBLDR — JOB 


MEMORY 


'I INCLUDING IOCS 

jnCK ORIGIN 

I 

jlBBF.R of' FILES - 

1 S.FBIN 
\ S.FBCU 

i PROGRAM 


12220 

12243 


■1 

DECK 

f 



12266 

1 

oecK 

*MIST 



35051 

4 ■ 

deck 

•HIX 



56135 

j 

'it * 

DECK 

' AMT 



56342 

i 

DECK 

'KKM : 



56546 

i 

deck 

ALE-C ■ 



63161 

1. 

oer.K 

’TICK 



63704 

1.. 

DECK 

* FUN SON 



64371 

I.' 

DECK 

»SST '■ 

f . 


65025 

I.-' 

SUBR 

» insyfb 



65140 

i. 

SURR 

•OUSYFS 


■ : 1 

65177 

h'- 

SUBR 

•POSTX 



65230 

- k'- 

SUBR 

•CNSTNT 



65542 

SUBR 

,• FPR 



65551 


SUBR 

• FRC 



65552 


SURF. 

•ICS 



65553 

jj,' 

SUBR 

•RWD 

t , 


66032 


SUBR 

•FXV 

f 


67206 

1* 

SUBR 

• FCV 



674 54 

SUBR 

•HCV 



67546 

li ■ 

SUBR 

•ICV 



67651 

a. 

SUBR 

•xcv' 

t ■ 


67671 

th 

SUBR 

• INT J 



67707 

' t4. 

SUBR 

• FFC 



70223 


SUBR 

*fpt 



70645 


00000 THRU 12211 
12220 


12266 THRU 71713 



